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Abstract  of  Dissertation  Presented  to  the  Graduate  Council 
of  the  University  of  Florida  in  Partial  Fulfillment  of  the  Requirements 
for  the  Degree  of  Doctor  of  Philosophy 


THE  MOBILITY,  RESISTIVITY  AND  CARRIER  DENSITY 
IN  p-TYPE  SILICON  DOPED  WITH  BORON,  GALLIUM  AND  INDIUM 


Luis  Carlos  Linares 
August  1979 

Chairman:  Sheng-San  Li 

Major  Department:  Electrical  Engineering 

Using  the  relaxation  time  approximation  and  a  three-band  model 
(i.e.,  nonparabolic  light-hole  band,  parabolic  heavy-hole  and  split-off 
bands),  a  derivation  involving  the  use  of  the  Boltzmann  transport  theory 
was  applied  to  obtain  expressions  for  the  valence  band  density-of-states 
effective  mass,  m^,  the  valence  band  conductivity  effective  mass,  m£, 
and  the  valence  band  Hall  effective  mass,  mjlj,  of  holes  in  p-type  sili¬ 
con,  Values  of  effective  mass  calculated  from  this  model  reveal  the 
temperature  and  dopant  density  dependence  of  the  effective  mass  due  to 
the  nonparabolic  shape  of  the  bands.  With  these  values  of  effective 
mass  and  the  three-band  model,  theoretical  calculations  of  hole  mobil¬ 
ity,  resistivity,  Hall  factor  and  Hall  mobility  were  conducted  for 
p-type  silicon  doped  with  boron,  gallium  and  indium  for  dopant  densities 
from  1014  to  1018  cm”8  and  temperatures  between  100  and  A00  K.  Scatter¬ 
ing  contributions  from  acoustical  and  optical  phonons,  and  ionized  and 
neutral  impurities  were  considered.  In  addition,  interband  scattering 


t  A 
JL  U 


158 


for  the  case  of  acoustical  phonons,  the  effect  of  hole-hole  scattering 

on  the  various  scattering  mechanisms,  and  the  nonparabolicity  of  the 

valence  band  were  also  taken  into  account  in  the  calculations.  The 

valence  band  density-of-states  effective  mass,  mg,  was  found  to  vary 

from  0.6567  mQ  at  100  K  to  0.8265  mQ  at  400  K,  while  the  valence  band 

conductivity  effective  mass,  mg,  increased  from  a  value  of  0.3604  mQ 

at  100  K  to  a  value  of  0.4910  nr  at  400  K.  The  valence  band  Hall  effec- 

o 

tive  mass,  m£,  varies  from  0.2850  nr  at  100  K  to  0.5273  in  at  400  K. 

The  masses  mg  and  mg  showed  little  change  with  dopant  density,  but 
mg  varied  by  as  much  as  63  percent  at  100  K  over  the  range  of  dopant 
densities  considered.  The  Hall  factor  was  expressed  in  terms  of  a  mass 
anisotropy  factor  and  a  scattering  factor.  These  two  components  of  the 
Hall  factor  were  separately  evaluated  to  emphasize  their  individual 
contributions  to  the  Hall  factor.  Theoretical  values  of  the  Hall  factor 
vary  between  1.73  and  0.77  over  the  dopant  density  and  temperature 
ranges. 

To  verify  our  theoretical  calculations,  resistivity  and  Hall 

coefficient  measurements  were  performed  on  silicon  samples  doped  with 

15 

boron,  gallium  and  indium  with  dopant  densities  ranging  from  4.25x10 

17 

to  9.07x10  cm  ,  for  100  s  T  s  400  K,  using  planar  square  array  test 

structures.  Agreement  between  the  calculated  and  measured  resistivity 
values  was  within  10  percent  over  the  range  of  temperatures  and  dopant 
density  studied.  Agreement  between  our  calculated  and' experimental 
values  of  Hall  mobility  was  within  15  percent  for  temperatures  greater 
than  150  K.  At  300  K,  agreement  between  theoretical  values  of  Hall 
factor  and  chose  deduced  from  experimental  data  of  this  work  and  that 
of  other  studies  was  within  15  percent  for  dopant  densities  greater 
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than  5x10*®  cm"’*.  From  the  results  of  this  study,  we  conclude  that  the 
theoretical  model  and  expressions  developed  here  can  accurately  describe 
the  mobility  and  resistivity  in  p-type  silicon.  The  omission  of  band 
anisotropy  considerations,  however,  introduces  considerable  error  in  the 
evaluation  of  the  Hall  factor  for  low  dopant  densities. 
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Using  the  relaxation  time  approximation  and  a  three-band  model 
(i.e.,  nonparabolic  light-hole  band,  parabolic  heavy-hole  and  split-off 
bands),  a  derivation  involving  the  use  of  the  Boltzmann  transport  theory 
was  applied  to  obtain  expressions  for  the  valence  band  density-of-states 
effective  mass,  m*,  the  valence  band  conductivity  effective  mass,  m^, 
and  the  valence  band  Hall  effective  mass,  mpj,  of  holes  in  p-type  sili¬ 
con.  Values  of  effective  mass  calculated  from  this  model  reveal  the 
temperature  and  dopant  density  dependence  of  the  effective  mass  due  to 
the  nonparabolic  shape  of  the  bands.  With  these  values  of  effective 
mass  and  the  three-band  model,  theoretical  calculations  of  hole  mobil¬ 
ity,  resistivity,  Hall  factor  and  Hall  mobility  were  conducted  for 
p-type  silicon  doped  with  boron,  gallium  and  indium  for  dopant  densities 
from  1014  to  10'®  cm’^  and  temperatures  between  100  and  400  K.  Scatter¬ 
ing  contributions  from  acoustical  and  optical  phonons,  and  ionized  and 
neutral  impurities  were  considered.  In  addition,  interband  scattering 
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for  the  case  of  acoustical  phonons,  the  effect  of  hole-hole  scattering 

on  the  various  scattering  mechanisms,  and  the  nonparabolicity  of  the 

valence  band  were  also  taken  into  account  in  the  calculations.  The 

valence  band  density-of-states  effective  mass,  mg,  was  found  to  vary 

from  0.6567  nv  at  100  K  to  0.8265  in  at  400  K,  while  the  valence  band 
o  o 

conductivity  effective  mass,  mg,  increased  from  a  value  of  0.3604  mQ 
at  100  K  to  a  value  of  0.4910  mQ  at  400  K.  The  valence  band  Hall  effec¬ 
tive  mass,  m£,  varies  from  0.2850  in  at  100  K  to  0.5273  itl  at  400  K. 

The  masses  mg  and  mg  showed  little  change  with  dopant  density,  but 
mg  varied  by  as  much  as  63  percent  at  100  K  over  the  range  of  dopant 
densities  considered.  The  Hall  factor  was  expressed  in  terms  of  a  mass 
anisotropy  factor  and  a  scattering  factor.  These  two  components  of  the 
Hall  factor  were  separately  evaluated  to  emphasize  their  individual 
contributions  to  the  Hall  factor.  Theoretical  values  of  the  Hall  factor 
vary  between  1.73  and  0.77  over  the  dopant  density  and  temperature 
ranges. 

To  verify  our  theoretical  calculations,  resistivity  and  Hall 

coefficient  measurements  were  performed  on  silicon  samples  doped  with 

boron,  gallium  and  indium  with  dopant  densities  ranging  from  4.25x1 0 ' ^ 
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to  9.07x10  cm  ,  for  100  <;  T  s  400  K,  using  planar  square  array  test 

structures.  Agreement  between  the  calculated  and  measured  resistivity 
values  was  within  10  percent  over  the  range  of  temperatures  and  dopant 
density  studied.  Agreement  between  our  calculated  and'  experimental 
values  of  Hall  mobility  was  within  15  percent  for  temperatures  greater 
than  150  K.  At  300  K,  agreement  between  theoretical  values  of  Hall 
factor  and  those  deduced  from  experimental  data  of  this  work  and  that 
.  of  other  studies  was  within  15  percent  for  dopant  densities  greater 
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than  5x10  cm  .  From  the  results  of  this  study,  we  conclude  that  the 
theoretical  model  and  expressions  developed  here  can  accurately  describe 
the  mobility  and  resistivity  in  p-type  silicon.  The  omission  of  band 
anisotropy  considerations,  however,  introduces  considerable  error  in  the 
evaluation  of  the  Hall  factor  for  low  dopant  densities. 
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Using  the  relaxation  time  approximation  and  a  three-band  model 
(i.e.,  nonparabolic  light-hole  band,  parabolic  heavy-hole  and  split-off 
bands),  a  derivation  involving  the  use  of  the  Boltzmann  transport  theory 
was  applied  to  obtain  expressions  for  the  valence  band  density-of-states 
effective  mass,  m*,  the  valence  band  conductivity  effective  mass,  m£, 
and  the  valence  band  Hall  effective  mass,  m^,  of  holes  in  p-type  sili¬ 
con.  Values  of  effective  mass  calculated  from  this  model  reveal  the 
temperature  and  dopant  density  dependence  of  the  effective  mass  due  to 
the  nonparabolic  shape  of  the  bands.  With  these  values  of  effective 
mass  and  the  three-band  model,  theoretical  calculations  of  hole  mobil¬ 
ity,  resistivity,  Hall  factor  and  Hall  mobility  were  conducted  for 
p-type  silicon  doped  with  boron,  gallium  and  indium  for  dopant  densities 
from  1014  to  10^®  cm“^  and  temperatures  between  100  and  400  K.  Scatter¬ 
ing  contributions  from  acoustical  and  optical  phonons,  and  ionized  and 
neutral  impurities  were  considered.  In  addition,  interband  scattering 
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for  the  case  of  acoustical  phonons,  the  effect  of  hole-hole  scattering 
on  the  various  scattering  mechanisms,  and  the  nonparabolicity  of  the 
valence  band  were  also  taken  into  account  in  the  calculations.  The 
valence  band  density-of-states  effective  mass,  mg,  was  found  to  vary 
from  0.6567  mQ  at  100  K  to  0.8265  mQ  at  400  K,  while  the  valence  band 
conductivity  effective  mass,  mg,  increased  from  a  value  of  0.3604  mQ 
at  100  K  to  a  value  of  0.4910  mQ  at  400  K.  The  valence  band  Hall  effec¬ 
tive  mass,  mj,  varies  from  0.2850  m„  at  100  K  to  0.5273  in  at  400  K. 

The  masses  mg  and  mg  showed  little  change  with  dopant  density,  but 
mg  varied  by  as  much  as  63  percent  at  100  K  over  the  range  of  dopant 
densities  considered.  The  Hall  factor  was  expressed  in  terms  of  a  mass 
anisotropy  factor  and  a  scattering  factor.  These  two  components  of  the 
Hall  factor  were  separately  evaluated  to  emphasize  their  individual 
contributions  to  the  Hall  factor.  Theoretical  values  of  the  Hall  factor 
vary  between  1.73  and  0.77  over  the  dopant  density  and  temperature 
ranges. 

To  verify  our  theoretical  calculations,  resistivity  and  Hall 

coefficient  measurements  were  performed  on  silicon  samples  doped  with 
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boron,  gallium  and  indium  with  dopant  densities  ranging  from  4.25x10 

17 

to  9.07x10  cm  ,  for  100  s  T  s  400  K,  using  planar  square  array  test 

structures.  Agreement  between  the  calculated  and  measured  resistivity 
values  was  within  10  percent  over  the  range  of  temperatures  and  dopant 
density  studied.  Agreement  between  our  calculated  and' experimental 
values  of  Hall  mobility  was  within  15  percent  for  temperatures  greater 
than  150  K.  At  300  K,  agreement  between  theoretical  values  of  Hall 
factor  and  those  deduced  from  experimental  data  of  this  work  and  that 
of  other  studies  was  within  15  percent  for  dopant  densities  greater 
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than  5x10  cm  .  From  the  results  of  this  study,  we  conclude  that  the 
theoretical  model  and  expressions  developed  here  can  accurately  describe 
the  mobility  and  resistivity  in  p-type  silicon.  The  omission  of  band 
anisotropy  considerations,  however,  introduces  considerable  error  in  the 
evaluation  of  the  Hall  factor  for  low  dopant  densities. 


CHAPTER  I 


INTRODUCTION 

The  goal  of  this  study  has  been  to  measure  and  compare  with  theory 
the  resistivity  and  Hall  mobility  of  holes  in  silicon  doped  with  gallium 
and  indium  as  functions  of  temperature  and  dopant  density.  Data  taken 
on  boron-doped  silicon  was  also  included  in  order  to  further  confirm  the 
adequacy  of  the  theoretical  model.  In  order  to  determine  theoretically 
the  resistivity  and  Hall  mobility,  one  must  first  calculate  the  conduc¬ 
tivity  mobility,  the  density  of  holes,  and  the  Hall  factor.  This  can 
only  be  done  with  a  thorough  understanding  of  the  energy  band  structure, 
the  scattering  mechanisms  involved,  and  the  carrier  statistics. 

The  application  of  an  electric  or  a  magnetic  field  to  a  crystal 
results  in  a  variety  of  carrier  transport  phenomena.  These  phenomena 
are  associated  with  the  motion  of  c  rrent  carriers  in  the  conduction  or 
valence  bands.  The  free  charge  carriers  in  a  semiconductor  will  acquire 
a  drift  velocity  under  the  influence  of  an  applied  electric  field.  This 
velocity  is  tne  net  result  of  the  momentum  gained  from  the  externally 
applied  field,  and  the  momentum  lost  in  collisions  which  tend  to  randomize 
the  carrier  momentum  [1].  If  the  field  is  expressed  in  volts  per  centi- 
meter,  and  the  velocity  in  centimeters  per  second,  a  mobility  is  defined 
as  the  incremental  average  speed  per  unit  electric  field,  and  is  expressed 
in  squared  centimeters  per  volt  second.  The  velocity,  and  consequently 
the  mobility,  is  determined  by  the  different  types  of  collisions  which 
the  carriers  undergo.  Collisions  of  carriers  with  lattice  atoms  which 
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are  out  of  their  equilibrium  positions  because  of  thermal  vibration, 
provide  an  upper  limit  to  the  mobility.  Scattering  of  the  carriers  is 
also  caused  by  impurities,  both  ionized  and  neutral.  At  high  tempera¬ 
tures  and  low  dopant  densities,  scattering  by  lattice  phonons  is  more 
effective  while  at  low  temperatures  and  high  impurity  densities,  scat¬ 
tering  by  ionized  and  neutral  impurities  predominates.  In  addition, 
the  effects  of  hole-hole  scattering  on  the  lattice  and  ionized  impurity 
scattering  mechanisms  need  to  be  considered.  Thus  in  calculating  the 
mobility  over  a  wide  range  of  temperatures  and  dopant  densities,  all  the 
different  scattering  mechanisms  must  be  taken  into  account.  Besides  the 
mobility,  the  density  of  holes  enters  the  problem  of  determining  the 
electrical  resistivity.  The  density  of  holes  associated  with  dopant 
atoms  is  a  function  of  the  ionization  energy  of  the  dopant  atom,  the 
temperature,  and  the  degeneracy  factor.  Calculation  of  the  Hall  mobility 
requires  knowledge  of  the  Hall  factor  which  is  a  function  of  the  scat¬ 
tering  mechanisms  and  effective  masses. 

For  purposes  of  device  design  it  is  necessary  to  know  the  correct 
relationship  between  the  resistivity  and  dopant  density  at  different 
temperatures.  Evaluation  of  the  characteristics  of  semiconductor  devices 
and  the  study  of  transport  phenomena  in  semiconductors  requires  an 
accurate  knowledge  of  variations  in  the  effective  mass,  mobility,  and 
carrier  density  with  changes  in  temperature  and  dopant  density.  Because 
of  this,  numerous  studies  of  mobility,  resistivity,  and  hole  density  in 

p-type  silicon  have  been  conducted  [1-16].  However,  due  to  the  complexity 
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of  the  valence  band  of  silicon  and  the  various  scattering  mechanisms 
involved,  these  studies,  for  the  most  part,  have  either  been  conducted 
in  temperature  and  dopant  density  ranges  designed  to  explore  only  a 
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particular  type  of  scattering  mechanism,  or  have  not  advanced  the  theory 

necessary  to  describe  the  experimental  result.  For  example,  Costato  and 

Reggiani  [4]  calculated  the  mobility  of  holes  for  pure  p-type  silicon  in 

which  lattice  scattering  dominates;  Braggins  [1]  considered  nonparaboli- 

city  and  all  the  relevant  scattering  mechanisms  with  the  exception  of 

hole-hole  scattering,  but  he  limited  his  investigation  to  dopant  densities 
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below  5x10  cm  and  low  temperatures;  Morin  and  Maita  [5]  considered 
wide  ranges  of  temperature  and  dopant  densities,  but  did  not  provide  a 
theoretical  examination  of  the  data.  Recently,  Li  [17]  developed  a 
theoretical  model  capable  of  describing  the  mobility  and  resistivity  of 
p-type  silicon  over  a  wide  range  of  temperatures  and  dopant  densities. 

This  improved  model  was  applied  to  the  case  of  boron-doped  silicon  with 
great  success  [17].  The  improvement  in  the  theory  consisted  mainly  of 
the  inclusion  of  hole-hole  scattering  effects,  and  consideration  of  the 
nonparabolic  nature  of  the  bands.  In  this  study,  Li's  model  [17]  has 
been  improved  by  including  consideration  of  interband  scattering  effects 
on  the  acoustic  phonon  scattering  mechanism,  and  has  been  applied  to  the 
study  of  silicon  doped  with  impurities  other  than  boron. 

With  some  exceptions  [14-16],  most  of  the  research  in  p-type  silicon 
has  been  conducted  with  boron  as  the  doping  impurity,  since  boron  is  the 
shallowest  acceptor  in  silicon  and  this  material  is  widely  available. 

A  very  limited  amount  of  data  is  available  on  silicon  doped  with  deeper 
impurities  such  as  gallium  and  indium.  These  dopants,  especially  indium, 
are  of  great  interest  to  modern  technology  because  of  their  application 
to  photo-detector  devices.  Curves  of  resistivity  and  mobility  as  func¬ 
tions  of  dopant  density  [2,3]  have  been  applied  to  characterizing 
boron-doped  starting  material  and  diffused  boron  layers  in  silicon,  and 
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have  been  found  highly  useful.  Similar  curves  developed  in  this  research 
may  be  expected  to  be  equally  useful  for  characterizing  and  integrating 
infrared  detectors  based  on  the  deeper  levels  of  indium  and  gallium  with 
on-chip  silicon  electronics.  Application  of  a  more  complete  theory  of 
mobility  and  resistivity  [17]  to  the  case  of  silicon  doped  with  gallium 
and  indium  should  provide  an  accurate  description  of  the  transport  of 
holes  in  this  material.  These  results  may  be  of  significant  use  in  the 
study  and  design  of  infrared  photo-detector  devices. 

In  this  research  the  mobility,  resistivity,  and  hole  density  have 

been  studied  over  a  temperature  range  from  100  to  400  K  and  dopant 

densities  from  4.25x10^  to  9.05x10^  cm"^.  Because  of  the  complexity 

brought  about  by  heavy  doping  effects  and  uncertainties  in  accounting 

for  hole  density  and  impurity  density  at  high  dopant  densities,  the 
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theoretical  analysis  has  been  restricted  to  densities  below  10  cm"  in 
which  the  use  of  Boltzmann  statistics  is  justified.  The  nonparabolic 
nature  of  the  valence  band  structure,  and  derivation  of  expressions  for 
the  temperature  dependent  effective  masses  are  presented  in  Chapter  II. 
Since  effective  mass  is  directly  related  to  the  shape  of  the  valence 
bands,  the  result  is  an  effective  mass  which  varies  with  temperature  and 
dopant  density.  The  mobility  formulation  includes  consideration  of  the 
relevant  scattering  mechanisms  and  how  these  are  modified  by  hole-hole 
scattering  effects.  These  scattering  mechanisms  are  considered  in  detail 
in  Chapter  III.  Since  the  different  scattering  mechanisms  which  contri¬ 
bute  to  the  mobility  have  different  temperature  and  energy  dependences, 
the  use  of  numerical  methods  and  curve  fitting  has  been  applied  in 
analyzing  the  data.  The  temperature  and  dopant  density  dependence  of 
resistivity  and  hole  density  is  analyzed  in  Chapter  IV.  In  Chapter  V, 


the  Hall  factor  is  discussed,  and  theoretical  calculations  of  Hall 
mobility  are  presented.  Fabrication  techniques  and  experimental  proce¬ 
dures  are  described  in  Chapter  VI.  Comparisons  of  experimental  results 
with  predictions  based  on  the  theory  of  Chapters  III  through  V  are  made 
in  Chapter  VII;  in  this  chapter  the  theoretical  results  are  also  compared 
with  data  published  by  other  workers.  Chapter  VIII  summarizes  the 
research  and  states  the  main  conclusions  derived  from  this  work. 


CHAPTER  II 


BAND  STRUCTURE  AND  EFFECTIVE  MASS 

2.1  Introduction 

The  interpretation  of  transport  properties  in  silicon  and  the  model¬ 
ing  of  silicon  junction  devices  depend  on  an  accurate  knowledge  of  values 
of  effective  mass.  The  complex  valence  band  structure  of  silicon  leads 
to  difficulties  in  the  study  of  transport  properties  of  holes  in  this 
material.  Thus  the  development  of  a  model  incorporating  the  nonparabolic 
nature  of  the  band  into  a  single  parameter,  the  combined  hole  effective 
mass,  would  greatly  simplify  the  study  of  mobility,  resistivity,  and  the 
Hall  effect  in  silicon.  Including  the  band  nonparabolicity  in  calcula¬ 
tions  of  relaxation  time  via  the  effective  mass  formulation  is  a  reason¬ 
able  procedure  and  has  been  applied  effectively  by  Radcliffe  [18]  to 
study  acoustic  phonon  scattering,  and  by  Barrie  [19]  to  study  optical 
phonon  and  impurity  scattering  in  nonparabolic  bands.  In  this  chapter 
we  will  derive  such  a  theoretical  model  for  hole  effective  mass  calcula¬ 
tions  in  silicon. 

Lax  and  Mavroides  [20]  have  derived  expressions  for  density-of-states 
effective  masses  and  mj^  for  the  heavy-hole  band  and  the  light-hole 
band,  respectively,  which  lead  to  the  generally  accepted  and  quoted 
value,  mg  =  0.591  mQ.  This  value,  however,  can  only  be  considered 
applicable  at  4.2  K,  where  m|  =  0.537  mQ  and  =  0.153  mQ.  A  number 
of  experimental  data  has  been  published  which  indicates  both  electron 
and  hole  effective  mass  to  be  dependent  both  on  temperature  and  dopant 
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density  [21,22].  Below  50  K,  hole  effective  mass  remains  constant  as 
indicated  in  high  frequency  magnetoconductivity  experiments  [23]. 

However,  at  higher  temperatures  and  for  higher  acceptor  impurity  densi¬ 
ties,  two  mechanisms  are  responsible  for  the  temperature  dependence  of 
the  effective  mass:  the  thermal  expansion  of  the  lattice,  and  the 
explicit  effect  of  temperature.  The  effect  of  the  thermal  expansion 
can  be  estimated  from  the  stress  dependence  of  the  effective  mass  [24], 
and  has  been  shown  to  be  negligible  [21,25].  The  explicit  temperature 
effect  however,  is  of  great  importance.  It  consists  of  three  parts: 

(a)  the  temperature  variation  of  the  Fermi  distribution  function  in  a 
nonparabolic  band,  (b)  the  temperature  dependent  distribution  function 
of  the  split-off  band,  and  (c)  the  temperature  variation  of  the  curva¬ 
ture  at  the  band  extremum  due  to  the  interaction  between  holes  and 
lattice  phonons. 

Following  the  work  of  Lax  and  Mavroides  [20],  but  using  Fermi-Dirac 
statistics  and  a  simplified  model  of  the  valence  band  structure  for 
silicon.  Barber  [25]  obtained  an  expression  for  the  density-of-states 
effective  mass,  mg,  which  is  temperature  and  hole-density  dependent. 
Barber,  however,  did  not  apply  the  nonparabolic  model  of  the  valence 
band  to  the  study  of  conductivity  or  Hall  effective  mass  in  p-type 
silicon.  Costato  and  Reggiani  [26]  also  developed  expressions  for  mg 
and  mg,  the  band  conductivity  effective  mass,  which  show  a  variation 
with  temperature,  but  they  neglected  the  effects  of  the  split-off  band 
and  the  temperature  variation  of  the  band  curvature. 

In  this  study,  the  expressions  for  density-of-states  effective  mass, 
conductivity  effective  mass,  and  Hall  effective  mass  of  holes  are  derived 
based  on  the  following  definitions.  The  density-of-states  effective 
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mass,  mg,  enters  in  the  normalization  of  the  distribution  function;  the 
conductivity  effective  mass,  mg,  is  the  mass  of  a  mobile  charge  carrier 
under  the  influence  of  an  external  electric  field;  and  the  Hall  effec¬ 
tive  mass,  mg,  is  the  mass  of  a  mobile  charge  carrier  under  the  applica¬ 
tion  of  external  electric  and  magnetic  fields.  The  reason  for  these 
particular  definitions  of  effective  masses  is  that  the  primary  applica¬ 
tion  of  this  work  is  to  generate  improved  theoretical  calculations  of 
Hall  mobility,  resistivity,  and  conductivity  mobility  [17].  The 
derived  expressions  were  used  to  calculate  hole  effective  masses  in 
p-type  silicon  over  a  wide  range  of  temperature  and  dopant  density. 

Since  the  crystal  structure  of  silicon  has  cubic  symmetry,  the  ohmic 
mobility  and  the  low-field  Hall  coefficient  are  isotropic.  An  angular 
average  of  the  effective  masses  may  be  performed  taking  into  account 
separately  the  warping  of  the  individual  bands  so  that  expressions  for 
mg,  mg,  and  m*  of  isotropic  form  can  be  derived.  Values  calculated  from 
these  expressions  differ  from  one  another  because  of  the  warping  and 
nonparabolicity,  and  consequently  effective  mass  in  each  band  depends  on 
temperature  and  dopant  density  in  its  own  way.  The  valence  band  struc¬ 
ture  of  silicon  is  presented  in  Section  2.2,  and  in  Section  2.3  expres¬ 
sions  for  mg,  mg,  and  m*  are  derived. 

2.2  The  Valence  Band  Structure  of  Silicon 

Theoretical  calculations  by  Kane  [27]  have  established  some  basic 
features  of  the  valence  band  of  silicon.  It  consists  of  heavy-hole  and 
light-hole  bands,  degenerate  at  k  =  0,  and  a  third  band  displaced  down 
in  energy  at  k  =  0  by  spin  orbit  coupling. 

The  heavy-hole  band  is  characterized  by  holes  with  an  energy  inde¬ 
pendent,  but  direction-dependent  effective  mass.  The  light-hole  band 
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is  characterized  by  holes  with  an  energy  and  direction-dependent  effec¬ 
tive  mass.  These  two  bands  can  be  described  by  the  E  vs  k  relationship 
[28] 


E(k) 


(2.1) 


where  A,  B,  and  C  are  the  experimentally  determined  inverse  mass  band 
2  2  2  % 

parameters,  k  =  (kx  +  k^  +  k2)  ,  and  the  upper  sign  is  associated  with 
the  holes  in  the  light-hole  band,  while  the  lower  sign  is  associated 
with  the  holes  in  the  heavy-hole  band.  Values  of  A,  B,  and  C  are 
obtained  by  cyclotron  resonance  measurements  at  4  K  [22,29]. 

Although  warped,  the  bands  are  parabolic  for  small  values  of  IT. 
However,  for  larger  values  of  k,  the  bands  become  nonparabolic,  and  along 
the  <100>  and  <1 1 1 >  directions  the  heavy-  and  light-hole  bands  are 
parallel  over  most  of  the  Brillouin  zone.  This  situation,  however,  is 
not  strictly  valid  for  general  directions  [30].  The  assumption  of 
overall  parallelism,  while  questionable  in  III-V  compounds,  is  reasonable 
in  the  case  of  Ge  and  Si  [27,31].  The  split-off  band  is  separated  at 
k  =  0  by  an  energy  A  =  0.044  eV.  [32],  and  is  characterized  by  an  effec¬ 
tive  mass  which  is  independent  of  energy  and  direction.  If  the 
anisotropy  is  small,  the  square  root  in  equation  (2.1)  may  be  expanded 
[20]  and  the  energy  surfaces  may  be  expressed  by 


Vr-  (A  ±  B'  )j  (0,4>) 


(2.2) 


where 
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B'  =  (B2  +  C2/6)?i  (2.3) 

8  and  <p  are  the  spherical  coordinates,  Ey  is  the  top  of  the  valence 
band,  and 


j(0,4>)  =  1  +  Y[sin4Q(cos4<i>  +  sin%)  +  cos48  -  2/3]  (2.4) 


with 


Y  =  ?  C2/2B ' (A  ±  B‘ )  (2.5) 

Following  the  work  of  Barber  [25],  we  have  used  the  simplified 
model  of  the  band  structure  illustrated  in  Figure  2.1.  In  this  model 
the  heavy-hole  band  is  considered  parabolic  and  thus  the  mass  m|  is  a 
constant,  equal  to  its  value  at  4.2  K.  For  energies  within  0.02  eV 
the  light-hole  band  is  considered  parabolic  with  a  constant  slope 
corresponding  to  the  value  of  m£  at  4.2  K.  For  higher  energies  the 
light-hole  band  is  assumed  to  take  on  approximately  the  same  slope  as 
that  of  the  heavy-hole  band,  but  remains  separated  from  the  heavy-hole 
band  by  A/3  eV  [27].  The  extrapolation  of  these  two  constant  slopes 
creates  the  kink  in  the  light-hole  band  at  0.02  eV.  Because  of  the 
change  in  slope,  the  light-hole  band  has  an  energy-varying  effective 
mass  and  in  general  can  only  be  described  in  terms  of  partial  Fermi- 
Dirac  integrals  [25].  Although  the  split-off  band  is  parabol ical ly 
distributed,  the  apparent  effective  mass  at  the  top  of  the  valence  band 
is  a  function  of  temperature  due  to  the  energy  displacement  at  k  =  0. 
Theoretical  and  exoerimental  studies  [33,34]  have  shown  that  at  high 
temperatures  the  heavy-hole' band  is  not  parabolic  and  thus  m|  is  not 
energy  and  temperature  independent.  However,  within  the  range  of 


Energy  (eV) 


Figure 


2.1.  Simplified  valence  band  structure  of  silicon  based  on 
Kane’s  [27]  calculations  and  measured  properties  of  the 
valence  band. 
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temperatures  considered  here,  the  assumption  of  parabolicity  for  the 
heavy-hole  band  based  on  Kane's  model  [27]  is  reasonable.  Other  studies 
[35,36]  support  the  validity  of  this  model  for  the  valence  band  of 
silicon. 


2.3  Effective  Mass  Formulation 

In  the  case  of  spherically  symmetric  energy  surfaces  all  of  the 
carriers  respond  in  the  same  way  to  a  given  set  of  applied  forces.  The 
effective  mass  then  acts  as  a  scalar  and  thus  has  the  same  value  for 
the  Hall  effect,  conductivity,  and  density  of  states.  For  nonspherical 
energy  surfaces,  however,  this  is  not  the  case.  The  mixed  response  of 
carriers  to  a  set  of  applied  forces  is  reflected  in  differences  between 
the  different  kinds  of  effective  masses.  The  densi ty-of-states  effective 
mass,  m*.,  is  defined  from  the  relationship 


p  =  — 

1  /F 


2TTkoTmDi 


3/2 


Fl/2(n) 


(2.6) 


where 


F 


1/2 


(n)  =  f 
o 


e'^de 

1  +  exp(c-n) 


(2.7) 


e  =  ( Ey  -  E)/kQT,  n  =  (Ey  -  Ep)/kQT,  kQ  is  the  Boltzmann  constant,  Ey  is 
the  top  of  the  valence  band,  and  i  =  1,  2,  3  refers  to  the  heavy-hole, 
light-hole,  and  split-off  bands,  respectively. 

The  electric  current  density  in  the  presence  of  electric  and  magnetic 
fields  can  be  expressed  by  [20] 

Jj  3  °jk£k  +  °jk£EkH£  +  ajk*i/k  Vm  + 


(2.8) 
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where  H^,  Hf()  are  the  electrical  and  magnetic  field  components  and 
the  o's  represent  single-energy-surface  conductivity  coefficients.  The 
first  coefficient  in  equation  (2.8)  is  the  zero-magnetic  field  electri¬ 
cal  conductivity,  and  the  second  coefficient  is  associated  with  the 
nondirectional  Hall  effect.  In  the  limit  of  weak  fields  the  expansion 
can  be  limited  to  the  first  two  terms.  We  use  the  electrical  conduc¬ 
tivity  coefficient,  o^.,  to  define  the  conductivity  effective  mass  m^, 
by  the  relationship 


Ci 


=  o 


Hi 


=  Pi 


e2<xi> 

mCi 


(2.9) 


and  the  Hall  effect  coefficient,  o^,  to  define  the  Hall  mobility  effec 
tive  mass  by  means  of  [37] 


°Hi  =  °123i  =  pi 


3  2 
eJ<tS 


(2.10) 


To  solve  for  mg^ ,  mg.,  and  mg.,  equations  (2.6),  (2.9)  and  (2.10)  are 
equated  to  the  following  expressions  for  p^,  and  o^: 


(2.11) 


4A2 


/  T 


8f0  D E  3 E  .3, 
3E  3kj  3kk  Q 


(2.12) 


rjk4 


4.V 


,  3fo  3E  3E  3 

i  T  ir  sir  nr  3ir 


{-  iL-i 
3kJ 


'£pq 


d3k 


(2.13) 
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where  f  is  the  Fermi-L)irac  distribution  function  and  c.  is  the  permu- 
o  Q 

tation  tensor.  Since  equations  (2.11)  through  (2.13)  do  not  assume  an 
effective  mass,  they  are  valid  both  for  parabolic  and  nonparabolic  band 
structures.  These  equations  are  then  evaluated  for  the  model  described 
in  Section  2.2. 

This  procedure  yields  single  m*^ ,  m^ ,  and  mjjj.  for  an  equivalent 
model  which  is  isotropic  and  parabolic.  These  values,  in  general,  will 
be  temperature  and  carrier-concentration  dependent.  Although  equations 
(2.6)  and  (2.11)  through  (2.13)  are  expressed  in  terms  of  Fermi-Dirac 
statistics  to  stress  their  generality,  conductivity  and  Hall  effective 
masses  were  derived  using  Boltzmann  statistics  to  simplify  the  form  of 
the  equation.  To  obtain  values  of  m^  and  m*^  we  also  require  a  proce- 
dure  for  evaluating  <x^>  and  <t^>  in  equations  (2.9)  and  (2.10).  This 
will  be  discussed  in  Chapter  III.  The  following  sections  present  the 
expressions  for  the  effective  masses  in  the  individual  bands. 


2.3.1  The  Heavy-Hole  Band 


In  this  band,  the  effective  masses  are  given  by 


7  mo  v  [f(-v)l2/3 

Ta^T 


(2.14) 


(2.15) 

(2.16) 


where  y  is  defined  in  equation  (2.5).  In  these  equations 
f (y)  =  (1  +  0.05y  +  0.01635y2  +  0.000908y3  +  ...) 


lb 


f-j  (y)  =  (1  +  0.01667y  +  0.041369y2  +  0.00090679y3 
+  0.00091 959y4  +  ...) 


and 

f2(y)  =  (1  -  0.01 667y  +  0.01 7956y2  -  C.0069857y3 
+  0.001 261 Oy4  +  ...) 

Since  the  heavy-hole  band  was  assumed  parabolic,  the  integrals  containing 
t  in  equations  (2.9)  and  (2.10)  are  identical  to  those  in  equations  (2.12) 
and  (2.13),  and  cancel  out. 


2.3.2  The  Liqnt-Hole  Band 

In  the  light-hole  band,  as  modeled  by  Figure  2.1,  the  effective 
masses  of  holes  are  obtained  in  terms  of  partial  Fermi -Dirac  integrals 
[25].  Thus 


(n^2) 


3/2  _ 


2m 


3/2 


/rT 


€AqT  , 

f(+T)  /  e'^de  + 

(A+B‘ )3/2  0  exp ( c ) 


f  (-Y  )n-j  00  cjde1 
(A-B 1 )3/2  l/kQJ  exp(el * 


(2.17) 
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3/2  . 

00  tot;  dc 

m?o  =  —r  I  ~r-nrs~  x 


2m. 


C2 


/F  o 


exp(e) 


f(+Y)  fS/koT  eV 
(A+B')3/2  o 


_ +  f(-y)ni  r  efdel 

(A-B')3/2  i/koT 


f^+Y)  f5/koT  T2c3/2de  f](-Y)n1  »  ^2el3/2dcl 

^  exp(e)  +  (A.b-)1/2  (/k  T  exp(£1 ) 


(A+B')1/2  o 


(2.18) 


(A+B ' )1/2f2(+Y)  /s/  0  x22e3/2exp(-e)dc  + 

o 


(A-B1 )]/2f2(-Y)n1  /  t A ,3/2exp(-e,)d£ 

2  '  C/k0T  2  ‘  1 


mH2  =  mo 


2_ 

/F 


f  ^3/2de  r  f (+Y )  _s2de 

o  exP(c)  L(A+B ' )3/2  o  eXp(£) 

f(-Y)nn 


e/k  t 


(A+B ' ) 

(A^B7?72  {/k  T  exP(c7r] 


Vt 


(2.19) 


where  e-j  =  e  -  A/3kQT,  c  -  C  -  A/3,  n-j  =  exp(-A/3kQT)  and  A  and  £  are 
defined  in  Figure  2.1 . 

In  this  case  because  equations  (2.11)  through  (2.13)  were 
expressed  in  terms  of  partial  Fermi-Dirac  integrals  and  equations  (2.6), 
(2.9)  and  (2.10)  were  expressed  in  terms  of  complete  Fermi-Dirac 
integrals,  the  dependence  on  x  does  not  cancel  out.  Thus  the  nonpara- 
bolicity  of  the  light-hole  band  introduces  a  dependence  on  the  scatter¬ 
ing  relaxation  time.  The  scattering  relaxation  time  is  discussed  in 
Chapter  III. 
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2.3.3  The  Split-Off  Band 

Although  the  split-off  band  is  parabolic,  the  apparent  effective 
mass  in  this  band  will  also  exhibit  a  temperature  dependence  due  to  the 
energy  displacement  at  k  =  0.  The  energy  of  a  hole  in  the  third  band 
is  given  by 

-2,2 

E  =  Ev-^-A-A  (2.20) 

o 


where  A  is  the  split-off  energy  (=  0.044eV),  and  A  is  one  of  the  inverse 
mass  band  parameters.  Substituting  equation  (2.20)  into  equations 
(2.11)  through  (2.13),  and  then  equating  to  equations  (2.6),  (2.9)  and 
(2.10)  for  the  split-off  band,  we  obtain 


m 


mD3  A°  exp  3k\} 


(2.21) 


ni  '  '3 

=  _o  o 


°o 

/  x^c3' 2exp(-e)dc 


C3  A  “  3/2 

J  x3e2-:>/  exp(-e2)de2 


(2.22) 


m 


|*  =  — h 

H3  A 


/  TT^c^//^exp(-e)dc 
o _ 

I  T3e23//2exp(-e2)de2 


(2.23) 


where  e2  =  c  -  A/kQT . 

The  combined  hole  densi ty-of-state  effective  mass  can  be  determined 
by  assuming  that  the  total  number  of  holes  in  the  valence  band  is  equal 
to  the  sum  of  the  holes  in  the  individual  bands 


t 
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p  =  p]  +  p2  +  p3 


(2.24) 


thus 


m*  =  C(m*1  )3/2  +  (m*2)3/2  +  (mfo)372]273  (2.25) 

This  combined  effective  mass  is  the  mass  corresponding  to  the  density- 

of-states  of  an  effective  single  equivalent  parabolic  valence  band. 

This  concept  is  useful  in  calculations  where  the  effective  density-of- 

states  at  different  temperatures  can  be  calculated  from  one  m*. 

The  explicit  temperature  variation  of  the  band  curvature  is  included 

by  assuming  that  the  density-of-states  near  the  band  edges  varies  in  a 

similar  manner  as  the  temperature  dependence  of  the  energy  gap  [25]. 

3/2 

Thus  (m^)  is  porportional  to  Eg  /Eg  where  EgQ  is  the  energy  gap  at  0  K. 

To  evaluate  the  total  band  equivalent  conductivity  and  Hall  effec¬ 
tive  masses,  we  assume  that  in  valence  band  conduction,  the  total  number 
of  holes  in  motion  is  equal  to  the  sum  of  the  holes  moving  on  the 
separate  energy  surfaces,  and  that  these  holes  can  be  modeled  as  moving 
on  a  single  spherical  energy  surface.  Thus,  the  ohmic  and  the  Hall 
conductivities  in  the  equivalent  valence  band  are  given  by 

aC  =  °C1  +  °C2  +  °C3  (2.26) 

and 

°H  =  °H1  +  °H2  +  °H3  (2.27) 


respectively. 
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Substituting  equations  (2.9)  and  (2.10)  into  equations  (2.26)  and 
(2.27)  it  follows  that 


A 

V 

rmBii 

1 3/2  i  "V 

[mD2l 

3/2  ,  <t3> 

'mD3j 

|3/2  1  1 

(  <T> 

imo  J 

1  mci  <T> 

*mD  > 

mC2  <T> 

-mD  - 

mC3  1 

(2.28) 

and 


m* 


r  <Ti  > 

'  l  <t2> 


3/2 
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1  ,  <t2  >  f  mD2 

2  2  !-3r- 
<T  > 


HI 


lmD 


3/2 


m; 


1  +<t32> 

T  + 


H2 


<x2> 


m; 


D3) 


D  J 


3/2 


m*  L 

mH3 


(2.29) 


Equations  (2.25),  (2.28)  and  (2.29)  were  evaluated  numerically  as  func¬ 
tions  of  temperature  and  acceptor  doping  density  for  p-type  silicon. 
Values  of  the  band  parameters,  j  A  j  =  4.27,  ) B  J  =  0.63  and  |C|  =  4.93, 
were  determined  at  4.2  K  by  Hensel  and  Feher  [22]  and  Balslev  and 
Lewaetz  [29].  In  order  to  simplify  the  calculations  and  maintain 
tractability,  anisotropies  in  the  relaxation  time  were  ignored.  A  rig¬ 
orous  analysis  of  the  conductivities  for  nonisotropic  scattering  would 
be  extremely  difficult  to  carry  out  because  no  relaxation  time  is 
expected  to  exist  in  the  usual  sense  [38]. 

Figure  2.2  shows  the  dependence  of  mj-j  with  temperature  in  the  range 
from  100  to  400  K.  The  slight  temperature  dependence  due  to  the  expli¬ 
cit  temperature  variation  of  the  curvature  at  the  edge  of  the  band 
results  in  an  effective  mass  increase  of  about  five  percent  in  each 
band  at  400  K,  This  can  be  seen  in  the  slope  of  mg-j .  The  temperature 
dependence  of  m^  is  more  pronounced  since  here  we  also  have  the 
effects  of  energy  displacement  at  k  =  0.  The  temperature  dependence  due 
to  nonparabolicity  is  very  apparent  in  the  shape  of  the  mj^  curve. 


Density-of-States  Effecti 
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The  temperature  dependence  of  the  conductivity  effective  mass  and 

the  Hall  effective  mass  is  shown  in  Figures  2.3  and  2.4,  with  the  dopant 
14  -3 

density  equal  to  10  cm  .  One  consequence  of  the  nonparabolicity  of 
the  light-hole  band  is  an  increase  in  the  valence  band  conductivity 
effective  mass  as  temperature  increases  from  100  to  400  K.  This  happens 
because  with  increasing  thermal  energy  kQT,  more  holes  reside  in  the  less 
parabolic  regions  of  the  light-hole  band.  The  results  plotted  in  Figure 
2.3  show  an  increase  in  m£  of  about  36  percent  in  this  temperature  range. 
The  temperature  dependence  of  m£j  can  be  attributed  mainly  to  the  non¬ 
parabolicity  of  the  light-hole  band.  In  the  temperature  range  from  100 
to  400  K,  m*  increases  from  0.2850  to  0.5273  m  .  The  slight  temperature 
dependence  of  m£-j  and  m*^  is  due  to  the  explicit  temperature  effect  and 
results  in  increases  of  7.7  percent  and  3.76  percent  in  the  m^-j  and  m*^ 
respectively.  A  larger  temperature  variation  occurs  in  the  case  of  the 
split-off  band  because  of  the  additional  effects  of  the  energy  displace¬ 
ment  at  k  =  0. 

Figures  2.5  and  2.6  show  the  variation  of  mjS  and  mj^  with  dopant 
density  and  temperature.  For  T  a  100  K,  m£  varies  less  than  10  percent 
in  the  dopant  density  range  from  10  ^  to  10^  cm”^.  Since  the  influence 
of  nonparabolicity  is  reduced  in  degenerate  material  [25],  it  follows 
as  shown  in  Figures  2.5  and  2.6  that  the  variation  of  effective  mass 
with  temperature  is  much  stronger  at  low  dopant  densities.  At  lower 
temperatures  there  is  a  much  greater  change  in  effective  mass  due  to 
variations  in  scattering  relaxation  time  with  percentage  of  ionized 
impurities. 


Temperature  (K) 


Figure  2.3.  Temperature  dependence  of  the  conductivity  effective  masses 
mCl  ’  mC2  and  mC3  ™  the  "'nd''vidua^  bands>  and  the  combined 

conductivity  effective  mass  m£  of  holes  in  silicon. 
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N.  =  10IH  cm 


acceptor  density  dependence  of  the  combined  conductivity  effective  mass 
ides  in  silicon  as  a  function  of  temperature. 


acceptor  density  dependence  of  the  combined  Hall  effective  mass  of  holes  in 
con  as  a  function  of  temperature. 
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2.4  Discussion 

The  idea  of  temperature-dependent  effective  mass  is  supported  by 
a  number  of  experimental  data.  Cardona  et  al  [21]  found  an  increase  of 
about  12  percent  in  optical  effective  mass  between  90  and  300  K  in 
heavily-doped  p-type  silicon.  Cyclotron-resonance  studies  conducted  by 
Hensel  and  Feher  [22]  show  that  when  carrier  heating  populates  deeper 
regions  of  the  light-hole  band,  the  nonparabolic  nature  of  this  band  at 
higher  values  of  k  results  in  an  increase  in  the  effective  mass  of  holes. 

The  model  used  here  in  the  calculation  of  hole  density-of-states 
effective  mass  is  identical  to  that  of  Barber  [25],  and  consequently 
our  results  for  mj*j.  and  m^  are  in  excellent  agreement  with  those  of 
Barber  [25].  We  have  extended  Barber's  work  to  the  calculations  of  mj~ 
and  mj^  in  p-type  silicon.  The  increase  of  m£  by  36  percent  at  400  K 
shown  in  Figure  2.3  is  much  larger  than  that  reported  by  Costato  and 
Reggiani  (9  percent)  [26].  Their  calculation  was  done  over  a  similar 
range  of  temperatures,  and  their  value  at  100  K,  mc=0.342  mQ,  is  some¬ 
what  lower  than  our  calculated  value  (.3604  mQ).  The  discrepancies 
between  our  results  and  those  of  Costato  and  Reggiani  are  due  mainly  to 
the  correction  of  for  the  explicit  temperature  dependence  of  the 
energy  gap,  the  inclusion  of  the  split-off  band,  and  the  consideration 
of  unequal  relaxation  times  in  the  three  bands.  Note  that  our  calcula¬ 
tions  of  effective  masses  were  achieved  through  more  rigorous  mathe¬ 
matical  derivations,  while  those  of  Costato  and  Reggiani  followed  a 
more  empirical  curve-fitting  type  of  procedure. 

The  experimental  values  of  density-of-states  effective  masses  of 
holes  in  p-type  silicon  have  been  published  by  numerous  authors  [21,22, 
39,40],  but  very  little  data  can  be  found  for  the  conductivity  and  the 
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Hall  effective  masses,  making  it  difficult  to  properly  assess  the  value 
of  our  calculations.  There  seems  to  be  no  obvious  way  to  measure  these 
quantities  from  d.c.  transport  measurements.  Magneto-kerr  effect 
measurements  conducted  by  Hauge  [41],  indicate  that  could  increase 
by  as  much  as  31  percent  in  the  range  of  temperatures  from  100  (M*  = 

0.510  Mq)  to  300  K.  This  is  in  reasonable  agreement  with  our  calculated 
percentage  increase  in  m*  in  the  same  temperature  range  (33  percent), 
but  it  is  impossible  to  compare  our  calculations  with  Hauge's  experi¬ 
mental  results,  because  our  effective  mass  definition  was  chosen  to  be 
mainly  applicable  to  the  study  of  the  Hall  and  conductivity  mobility  in 
the  low  field  limit,  and  this  may  not  apply  to  the  measurements  of 
Hauge  [41]. 

From  the  results  of  this  chapter  it  can  be  seen  that  the  approxima¬ 
tion  of  a  constant  effective  mass  seems  to  be  inadequate  to  describe 
transport  properties  of  holes  in  silicon  above  100  K.  There  is  a  sub¬ 
stantial  increase  in  the  effective  mass  of  holes  from  100  to  400  K  due 
to  the  nonparabolicity  of  the  light-hole  band,  and  a  smaller,  though  not 
negligible,  contribution  due  to  the  explicit  temperature  dependence  and 
the  effects  of  the  split-off  band.  The  validity  of  this  model  for  the 
calculation  of  densi ty-of-states  effective  mass  has  been  well  established 
[25].  Barber  [25]  has  shown  that  when  the  temperature-dependent  effec¬ 
tive  masses  are  substituted  into  the  theoretical  expression  for  intrinsic 
carrier  density  in  silicon,  the  agreement  with  reported  measurements  of 
n.j  is  within  the  limits  of  experimental  error.  Application  of  this 
model  to  theoretical  calculation  of  mobility  and  resistivity  in  p-type 
silicon  [17]  has  provided  excellent  agreement  between  theoretical  and 
and  experimental  values  (resistivity  with  ±6  percent)  over  a  temperature 
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range  from  100  to  400  K  and  dopant  density  range  from  10  to 
3x1018  cm-3.  This  calculation  is  limited  to  applications  in  conduc¬ 
tivity  mobility  and  low  field  Hall  effect. 


CHAPTER  III 


MOBILITY  AND  SCATTERING  RELAXATION  TIME 
3.1  Introduction 

The  study  of  transport  phenomena  in  semiconductors  requires  an 
accurate  knowledge  of  variations  in  the  conductivity  mobility  and  the 
resistivity  with  changes  in  temperature  and  dopant  density.  The  resis¬ 
tivity  is  an  easily-measured  parameter,  but  the  conductivity  mobility 
is  a  more  difficult  parameter  to  evaluate.  In  general,  four  different 
kinds  of  mobility  enter  into  common  discussion  [42].  The  microscopic 
mobility  is  the  actual  velocity  per  unit  electric  field  of  a  free 
carrier  in  a  crystal.  This  cannot  be  measured  directly.  The  conduc¬ 
tivity  mobility  is  the  mobility  associated  with  the  conductivity 
expression,  a  =  epp^,.  This  mobility  involves  an  average  relaxation 
time  <x>  dependent  on  the  nature  of  the  scattering  process,  and  in  the 
case  of  nonspherical  equal  energy  surfaces,  this  mobility  also  involves 
a  combined  effective  mass.  The  Hall  mobility  is  the  product  of  the 
measured  conductivity  and  the  measured  Hall  coefficient.  In  general, 
the  Hall  mobility  differs  from  the  conductivity  mobility  by  a  factor 
called  the  Hall  factor.  The  drift  mobility  is  the  velocity  or  drift 
per  unit  field  for  a  carrier  moving  in  an  electric  field.  If  trapping 
centers  are  present,  so  that  the  actual  drift  process  is  not  simply 
motion  through  the  conduction  band,  but  involves  a  series  of  trapping 
and  untrapping  processes,  the  drift  mobility  can  be  much  less  than  the 
conductivity  mobility.  The  four  mobilities  are  all  equal  only  when  the 
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following  three  conditions  are  met  [42]:  (a)  spherical  equal  energy 

surfaces  with  extremum  at  k  =  0,  (b)  relaxation  time  independent  of 
carrier  energy,  and  (c)  negligible  trapping  effects.  Since  conditions 
(a)  and  (b)  are  not  met  in  p-type  silicon,  it  is  improper  to  judge  the 
behaviour  of  one  kind  of  mobility  based  on  knowledge  of  a  different 
kind  of  mobility.  Thus  drift  or  Hall  mobility  data  cannot  be  tacitly 
assumed  to  be  accurate  substitutes  for  conductivity  mobility  values. 

As  mentioned  above,  the  conductivity  mobility  involves  an  average 
scattering  relaxation  time.  In  any  semiconductor,  the  charge  carriers 
(i.e.,  holes  and  electrons),  at  temperatures  above  absolute  zero,  may 
be  scattered  by  a  number  of  mechanisms.  Different  mechanisms  are 
dominant  in  certain  temperature  and  dopant  density  regimes,  but  in  some 
cases  two  or  more  may  be  interacting  simultaneously.  Thus  in  calculat¬ 
ing  the  conductivity  mobility  over  a  wide  range  of  temperatures  and 
dopant  densities,  all  the  relevant  scattering  mechanisms  must  be  taken 
into  account.  In  the  case  of  silicon,  acoustic  and  optical  phonon 
scattering,  and  ionized  and  neutral  impurity  scattering  are  of  major 
importance.  Hole-hole  scattering  also  plays  an  important  role  in  deter' 
mining  the  mobility.  In  the  following  sections  the  theoretical  effects 
of  these  scattering  mechanisms  on  the  mobility  will  be  considered. 

3.2  Mobility  and  Average  Scattering  Relaxation  Time 

The  calculation  of  mobility  of  holes  in  the  valence  band  of  silicon 
is  accomplished  by  evaluating  the  mobility  separately  in  the  heavy-hole 
band,  the  light-hole  band,  and  the  split-off  band  considering  all 
appropriate  scattering  mechanisms.  The  overall  mobility  is  then 
evaluated  as  a  weighted  average  of  the  single-band  mobilities  over  the 
individual  hole  densities  in  each  band. 


The  conductivity  mobility  in  each  of  the  three  valence  bands  is 
calculated  from 
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e<T^> 
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where 


for  the  case  of  Fermi-Dirac  statistics,  and  t-  represents  the  total 
scattering  relaxation  time  in  band  i.  Because  each  scattering  mechanism 
has  its  own  dependence  on  scattering  energy,  a  simple  closed  form 
expression  for  total  scattering  relaxation  time  as  a  function  of  temper¬ 
ature  cannot  be  obtained.  The  use  of  numerical  techniques  is  necessary 
to  solve  for  the  relaxation  time.  In  the  case  of  p-type  silicon,  the 
peculiarities  of  a  degenerate,  warped,  and  nonparabolic  valence  band 
must  be  taken  into  account  [1],  The  possibility  of  interband  as  well  as 
intraband  transitions  must  also  be  taken  into  account  in  the  analysis. 
With  the  inclusion  of  interband  scattering  as  given  by  Bir  et  al.  [43], 
the  total  relaxation  time  in  the  heavy-  (i  =  1)  and  light-holes  (i  =  2) 
bands  is  given  by 
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where 
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TH  “  [Taci”1  +  Vi"'  *  TIi"'  *  V"']  <3'5> 

The  total  relaxation  time  in  the  split-off  band  is  given  by 

x3  =  [xac3  1  +  xo3  1  +  TI3  1  +  XN3  ^  1  ^3-6^ 

Only  transitions  between  the  light-  and  heavy-hole  band  are  con¬ 
sidered;  the  relaxation  time  takes  into  account  a  transition  from 
band  i  to  band  j;  and  t^,  tqi->  and  are  the  relaxation  times 

corresponding  to  scattering  by  acoustical  phonons,  optical  phonons, 
ionized  impurities,  and  neutral  impurities  respectively,  with  i  as  the 
band  index.  The  procedure  for  including  the  nonparabol icity  of  the 
band  structure  into  calculations  of  relaxation  time,  consists  of 
modifying  the  relaxation  time  for  a  given  scattering  process  by  replac¬ 
ing  the  temperature  independent  effective  mass  of  the  parabolic  band  by 
the  temperature  dependent  effective  mass  of  the  nonparabolic  band. 

This  procedure  has  been  successfully  applied  to  the  study  of  acoustic 
phonon  scattering  in  nonparabolic  bands  by  Radcliffe  [18].  Optical 
phonon  and  ionized  impurity  scattering  in  nonparabolic  bands  have  been 
considered  by  Barrie  [19]  in  the  same  manner.  Braggins  [1]  has  used 
the  same  method  to  include  nonparabolicity  in  his  study  of  p-type  sili¬ 
con.  In  this  work,  the  relaxation  times  appropriate  to  degenerate, 
parabolic  valence  bands  have  been  used  and  modified  according  to  the 
prescription  of  Radcliffe  [18],  Barrie  [19],  and  Braggins  [1].  The 
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anisotropy  of  the  energy  spectrum  is  not  considered  in  this  model, 
because  from  the  transport  theory  for  parabolic  bands  it  is  known  that 
this  anisotropy  has  no  influence  on  the  temperature  dependence  of 
mobility,  but  only  on  its  absolute  value  [10].  Each  of  the  four 
scattering  mechanisms  will  now  be  discussed. 

3.3  Acoustical  Phonon  Scattering 

The  relaxation  time  for  scattering  by  acoustical  phonons  includes 
both  the  possibility  of  interband  as  well  as  intraband  scattering.  The 
treatment  of  the  acoustical  phonons  has  been  based  on  the  theory  of 
Bir,  Normantas,  and  Pikus  [43]  where  the  relaxation  times  can  be 
expressed  in  terms  of  a  single  constant,  x  ,  which  controls  the  overall 
magnitude  of  the  scattering.  Both  transverse  and  longitudinal  phonons 
participate  in  the  scattering  so  that 
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for  intraband  scattering,  while 


(3.9) 


for  interband  scattering.  In  the  split-off  band,  the  scattering  relaxa¬ 
tion  time  is  given  by 
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In  these  equations 
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Yij  *  mDi^m0j’  8  =  b/a»  a  a°d  b  are  valence  band  acoustic  deformation 
potential  constants  in  the  Picus  and  Bir  [44]  notation,  p$  is  the 
density,  and  Ct  are  the  longitudinal  and  transverse  sound  velocities 
in  silicon  and  L..  and  T. .  are  functions  of  6  and  Y- •  defined  in  [43]. 

*  J  '  J  I J 

3.4  Optical  Phonon  Scattering 

Optical  phonon  scattering,  while  negligible  at  very  low  tempera¬ 
tures,  cannot  be  ignored  at  high  temperatures.  Ehrenreich  and 
Overiiauser  [45]  have  calculated  the  mobility  of  holes  in  silicon  and  its 
dependence  on  temperature.  The  calculated  mobility  follows  a  T~2’^ 
dependence  for  reasonable  choices  of  the  parameters  which  described  the 
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mixing  of  optical  and  acoustical  phonon  scattering.  This  agrees  with 
experimental  results  [5,8].  The  relaxation  time  for  scattering  by 
nonpolar  optical  phonons  is  given  by  [46] 


where  is  the  Debye  temperature,  nQ  =  (exp(0p/T)-l )-^  is  the  phonon 
distribution  function,  and  W  is  a  constant  which  determines  the  rela¬ 
tive  coupling  strength  of  the  holes  to  the  optical  phonon  mode  compared 
to  the  acoustical  phonon  mode 
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where  Dq  is  the  optical  deformation  potential  constant.  The  first  term 
in  the  brackets  of  equation  (3.12)  corresponds  to  optical  phonon  emmis- 
sion  and  is  relevant  only  when  this  is  energetical ly  possible  (e>0p/T). 
The  second  term  in  the  brackets  corresponds  to  optical  phonon  absorption. 


3.5  Ionized  Impurity  Scattering 

The  Columbic  interaction  between  ionized  impurities  and  charge 
carriers  drifting  through  the  cyrstal  under  the  action  of  an  applied 
electric  field  causes  scattering  of  the  charge  carriers.  Scattering  by 
ionized  impurities  was  first  considered  by  Conwell  and  Weisskopf  [47]. 
The  basic  assumption  is  that  the  Coulomb  field  is  cut  off  at  half  the 
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distance  between  charged  impurities.  This  is  equivalent  to  assuming 
that  a  charge  carrier  sees  only  one  charged  impurity  at  a  time,  the 
effect  of  the  other  charged  impurities  being  sufficiently  screened  as 
to  be  negligible.  This  approach  was  improved  by  Brooks  [48]  and 
Herring  [49]  who  associated  the  cut-off  of  the  Coulomb  potential  with 
a  screening  distance,  the  free  carriers  being  assumed  to  provide 
screening  against  the  charge  of  the  impurities.  In  the  low  dopant 
density  limit,  the  scattering  relaxation  time  due  to  ionized  impurities 
is  given  by  [48,49] 
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where  p*  is  the  screening  carrier  density,  p'  =  p  +  NA(1  -  N^/N^),  for 

nd  *  0. 


3.6  Neutral  Impurity  Scattering 

Scattering  by  neutral  impurities  in  semiconductors  has  been  con¬ 
sidered  by  Erg insoy  [50]  as  a  variation  of  the  problem  of  the  scattering 
of  electrons  by  neutral  hydogen  atoms.  The  result  is  a  temperature 
independent  relaxation  time  given  by 
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where  is  the  density  of  neutral  impurities  and  is  the  geometric 
mean  mass  appropriate  for  evaluating  the  scaled  Bohr  radius  term  [48]. 
Sclar  [51,52]  has  included  the  possibility  of  bound  states  in  the 
evaluation  of  electron-hydrogen  impurity  scattering  by  using  a  three- 
dimensional  square  well  to  estimate  the  influence  of  a  weakly-bound 
state  on  the  scattering.  In  this  case  the  relaxation  time  is  given  by 
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is  the  binding  energy  of  neutral  acceptors. 

For  silicon  doped  with  shallow  impurities,  this  type  of  scattering 
is  important  at  low  temperatures  where  neutral  impurities  may  outnumber 
ionized  impurities.  For  the  deeper  levels,  where  neutral  impurities  can 
exist  at  higher  temperatures,  the  influence  of  neutral  impurity  scatter¬ 
ing  can  extend  over  a  wide  range  of  temperatures. 

3.7  Effect  of  Hole-Hole  Scattering 

The  expressions  thus  far  presented  for  scattering  relaxation  time 
neglect  the  effect  of  hole-hole  scattering.  Although  hole-hole  scatter¬ 
ing  does  not  affect  the  current  density  directly  since  it  cannot  alter 
the  total  momentum,  it  tends  to  randomnize  the  way  in  which  this  total 
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momentum  is  distributed  among  holes  of  different  energies.  When  the 
scattering  mechanism  is  such  as  to  lead  to  a  nonuniform  distribution,  * 
hole-hole  scattering  gives  rise  to  a  net  transfer  of  momentum  from 
holes  which  dissipate  momentum  less  efficiently  to  those  which  dissipate 
momentum  more  efficiently,  resulting  in  an  overall  greater  rate  of 
momentum  transfer,  and  lower  mobility  [53].  Thus  the  size  of  the  effect 
of  hole-hole  scattering  on  the  scattering  relaxation  time  is  a  function 
of  the  energy  dependence  of  the  relaxation  time.  The  hole-hole  reduc¬ 
tion  factor,  can  be  derived  by  means  of  a  classical  formulation 
introduced  by  Keyes  [54].  When  hole-hole  collisions  are  much  more 
frequent  than  hole-acceptor  collisions,  the  average  relaxation  time  for 
a  parabolic  band  in  the  Keyes  [54]  approximation  approaches  the 
limiting  form 


(3.20) 


where  f  is  the  Fermi-Dirac  distribution  function.  On  the  other  hand, 
o 

if  hole-hole  collisions  are  neglected,  the  average  relaxation  time  is 
given  by  equation  (3.2). 

Thus  the  hole-hole  reduction  factor  (i.e.,  the  ratio  of  <xhh>  to 
<t> )  can  be  expressed  as 
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for  optical  phonon  scattering,  and  y^  ,  the  hole-hole  reduction  factor, 
is  evaluated  from  equation  (3.21). 

For  acoustical  phonon  scattering  it  is  assumed  that  yhha  decreases 
linearly  with  increasing  dopant  density  from  a  value  of  one  to  a  value 
y^3  *  9n/32  «  0.88  [17]  in  a  certain  range  of  impurity  concentrati 
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The  exact  relationship  (y,  a  =  1.0004  -  4.013378  x  10'19NA,  1015  s  N.  s 
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3  s  10)  is  determined  empirically  with  a  best  fit  of  the  experimental 
data. 

Luong  and  Shaw  [55]  using  a  one-particle-like  approximation  from 
the  Hartree-Fock  theory,  have  shown  that  by  inclusion  of  hole-hole 
scattering,  the  Brooks-Herring  [48,49]  formula  is  reduced  by  a  factor 
which  can  be  expressed  in  closed  form  as 


(3.22) 


where  is  the  ionized  acceptor  density  and  p1  is  the  screening  hole 
density.  In  the  case  of  neutral  impurity  scattering,  hole-hole  scat¬ 
tering  has  no  significance  because  is  independent  of  hole  energy. 

Thus  the  overal  scattering  relaxation  time  in  each  hole  band  is 
calculated  from  equations  (3.3),  (3.5),  and  (3.6)  with  the  terms  of 
these  equations  properly. corrected  for  the  effects  of  hole-hole  scatter¬ 
ing.  Because  the  individual  energy  surfaces  are  different  from  each 
other,  the  relaxation  times  also  differ  from  each  other  and  cannot  be 
assumed  equal  except  in  restricted  ranges  of  temperature  and  dopant 
density  [43]. 
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3.8  Mobility  in  the  Combined  Valence  Band 

The  conductivity  mobility  in  each  individual  band  is  calculated 
from  equation  (3.1),  and  the  combined  conductivity  mobility  in  the 
valance  band  is  then  evaluated  as  a  weighted  average  of  the  single-band 
mobilities  over  the  population  of  holes  in  each  band,  thus 
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Using  equation  (3.23)  and  the  parameters  listed  in  Table  3-1,  we 

have  calculated  the  hole  mobility  for  silicon  doped  with  boron,  gallium, 

and  indium  as  functions  of  dopant  density  and  temperature,  for 

1014  s  <  1013cm"3  and  100  <  T  <  400  K.  The  results  are  displayed  in 

Figures  3.1  through  3.6.  In  the  calculations  of  mobility  and  resistivity 

in  silicon  doped  with  gallium  and  indium,  it  was  assumed  that  boron 

impurities  were  also  present.  Since  very  pure  silicon  has  a  resistivity 

1 3 

on  the  order  of  lOOQft-cm,  it  was  assumed  that  boron  densities  of  10 

13  -3 

and  5x10  cm  existed  in  the  gallium-  and  indium-doped  samples, 
respectively.  The  values  of  these  background  densities  were  deduced 
from  a  best  fit  of  the  experimental  data.  For  this  reason,  especially 
in  the  case  of  indium-doped  silicon,  the  actual  role  of  the  impurities 
at  low  temperatures  and/or  low  dopant  densities  is  masked  by  the  action 
of  the  always  present  boron  impurities.  As  the  dopant  density  and 
temperature  increase,  the  assumed  background  densities  of  boron 
impurities  become  insignificant  compared  to  the  density  of  ionized  dopant 
atoms,  and  Figures  3.1  through  3.6  accurately  depict  the  influence  of 
the  particular  type  of  impurity  on  the  resistivity  and  mobility  of  holes 
in  p-type  silicon.  The  figures  also  show  that  for  the  case  of  the 


Table  3-1.  Values  used  in  the  calculations. 


Parameter 

Value 

Unit 

A 

44.0 

meV 

a 

-6.4* 

eV 

b 

-1.36* 

eV 

C  2/C  2 
u z  /0t 

2.09* 

0D 

735 

K 

ps 

2.329  x  103 

kg/m3 

e„ 

11.7 

£ 

s 

0 

Tx 

6.96  x  10'10 

sec  K3/2 

W 

0.244 

mo 

9.1  x  10‘31 

kg 

h 

6.25  x  10'34 

joule-sec 

ko 

1.38  x  10'23 

joules/K 

e 

1.6  x  10"19 

coul 

*  These  values  were  obtained  from  references  [1]  and  [22]. 
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Figure  3.5.  The  calculated  hole  mobility  vs  temperature  for  gallium- 
doped  silicon  with  dopant  density  as  a  parameter. 
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Figure  3.6.  The  calculated  hole  mobility  vs  temperature  for 
doped  silicon  with  dopant  density  as  a  paramete 
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shallower  ionization  energies,  the  mobility  depends  more  strongly  on 
temperature  for  the  lightly-doped  case  where  lattice  scattering  is 
dominant  and  become  less  temperature  dependent  as  the  dopant  density 
increases. 

The  constant,  t  ,  was  found  by  fitting  the  mobility  to  experimental 
data  in  the  lattice-scattering-limited  range.  Our  value  of  x  is 
equivalent  to  an  acoustic  deformation  potential  constant  of  8.099  eV. 

The  optical  phonon  coupling  constant,  W,  was  then  found  by  fitting  the 
mobility  to  the  high  temperature  experimental  data.  Our  value  of  W  is 

g 

equivalent  to  an  optical  deformation  potential  constant  of  6.024x10  eV/ 
cm. 


CHAPTER  IV 


HOLE  DENSITY  AND  RESISTIVITY 


4.1  Introduction 

The  resistivity  of  semiconductor  materials  is  one  of  their  most 
useful  and  easily  measured  properties.  Theoretical  calculations  of 
resistivity  depend  on  the  formulation  of  conductivity  mobility,  and  the 
determination  of  hole  density.  For  extrinsic  semiconductors,  the  hole 
density  is  determined  primarily  by  the  percentage  of  ionization  of 
impurity  atoms.  The  following  sections  discuss  the  dependence  of  hole 
density  and  resistivity  on  temperature  and  dopant  density. 

4.2  Ionization  of  Impurity  Atoms 

For  the  case  of  Fermi-Dirac  statistics,  the  hole  density  is  given 
by 


where,  mg,  the  density-of-states  effective  mass,  contains  information 
pertaining  to  the  nonparabolic  nature  of  the  valence  band.  In  the  limit 
of  low  dopant  densities,  equation  (4.1)  reduces  to 

p  =  Nyexp(n)  (4.2) 

2  3/2 

where  Ny  =  2(2  tt  mg  kQT/h  )  is  the  effective  density  of  valence  band 
states.  For  the  range  of  temperatures  considered  in  this  study,  the 
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hole  density  is  calculated  by  assuming  that  the  density  of  carriers  is 
determined  by  the  impurities  present  in  the  silicon  sample.  The  density 
of  ionized  acceptor  impurities  in  p-type  silicon  is  computed  from  the 
charge  neutrality  equation 
-  + 

NA  -  ND  =  p  -  n  (4.3) 

This  reduces  to 

p  -  Na  (4.4) 


for  the  case  of  uncompensated  material. 

The  density  of  ionized  acceptors  is  [56] 


N 


A 


1  +  g  exp 


(ea-e 

ry 


(4.5) 


where  EA  is  the  acceptor  ionization  energy,  and  g  is  the  ground  state 
degeneracy.  Excited  states  have  a  very  minor  influence  on  the  carrier 
concentration  due  to  the  large  separation  between  the  ground  state  and 
the  excited  states  [1,56].  Letting 


g  =  4  +  2  exp 


(4.6) 


enables  us  to  include  the  contribution  of  the  split-off  band  [17].  The 
density  of  ionized  acceptors  is  computed  by  iterating  Ep  in  equations 
(4.2)  and  (4.5)  until  equation  (4.4)  is  fulfilled  within  a  given  level 
of  accuracy. 


Experimental  evidence  shows  that  the  acceptor  ionization  energy  E^ 

is  not  a  constant,  but  decreases  with  increasing  dopant  density  [9]. 

Penin  et  al.  [57]  have  determined  in  a  study  of  heavily  doped  silicon 

from  4  to  300  K  that  for  shallow  impurities  such  as  boron  and  phosphorus 

the  ionization  energy  decreases  and  finally  disappears  altogether  for 
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impurity  densities  greater  than  3x10  cm  .  For  impurities  with  deeper 
activation  energies,  it  is  also  expected  that  at  some  impurity  concen¬ 
tration,  the  impurity  activation  energy  should  become  a  function  of  the 
impurity  concentration.  However,  in  the  case  of  gallium  and  indium, 
this  should  happen  at  higher  impurity  concentrations  than  for  the 
shallower  level  impurities.  This  is  due  to  the  smaller  geometrical 
dimensions  of  the  wave  functions  applicable  to  the  deeper  levels,  so 
that  overlapping  effects  which  promote  the  reduction  in  activation  energy 
require  higher  impurity  concentrations  [14].  For  shallow  impurities 
such  as  boron  and  phosphorus,  empirical  expressions  [9,57]  relating  the 
dependence  of  ionization  energy  to  dopant  density  have  been  established. 
In  the  case  of  Ga,  there  is  data  [15]  on  activation  energy  vs  concentra¬ 
tion,  but  not  enough  on  which  to  base  an  accurate  relationship.  For 
this  reason  the  value  of  EA  =  0.056eV  was  used.  For  In,  =  0.156eV 
[58]  was  used.  Figures  4.1  through  4.3  show  the  ratio  of  ionized  and 
total  impurity  density  as  a  function  of  impurity  density  with  tempera¬ 
ture  as  a  parameter  for  100  s  T  s  400  K  for  silicon  doped  with  boron, 
gallium,  and  indium.  It  is  clearly  shown  in  these  figures  that  the 
ionization  of  impurities  for  the  deeper  levels  is  significantly  lower 
even  at  low  dopant  densities  so  that  it  is  necessary  to  go  to  higher 


I  <*«*-*- 


cal  calculations  of  the  ratio  of  ionized  and  total  indium  density 
im  density  with  temperature  as  a  parameter. 
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temperatures  to  achieve  total  ionization  of  impurities.  The  deioniza¬ 
tion  of  impurities  is  mos't  significant  for  low  temperatures  and  high 
impurity  densities. 


4.3  Resistivity  of  p-Type  Silicon 


The  resistivity  of  p-type  silicon  is  given  by 


(4.7) 


where  is  the  hole  conductivity  mobility  calculated  from  equation 
(3.23)  and  p  is  the  hole  density  discussed  in  Section  4.2.  Equation 
(4.7)  was  used  to  calculate  the  hole  resistivity  for  silicon  doped  with 
boron,  gallium,  and  indium  as  a  function  of  dopant  density  and  tempera¬ 
ture,  for  1014  <  <  1Q^®  cm”3  and  100  s  T  s  400  K.  The  results  are 

displayed  in  Figures  4.4  through  4.9.  In  the  calculations  of  resistiv¬ 
ity  in  silicon  doped  with  gallium  and  indium,  as  was  done  for  conductiv¬ 
ity  mobility,  it  was  assumed  that  boron  impurities  were  also  present. 
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Boron  densities  of  10  and  5x10  cm  were  assumed  to  exist  in  the 
gallium-  and  indium-doped  samples,  respectively.  The  values  of  these 
background  densities  were  deduced  from  a  best  fit  of  the  experimental 
data.  As  the  dopant  density  and  temperature  increase,  the  assumed 
background  densities  of  shallow  impurities  becomes  insignificant 
compared  to  the  density  of  ionized  dopant  atoms,  and  Figures  4.4  through 
4.9  accurately  depict  the  influence  of  the  particular  type  of  impurity 
on  the  resistivity  of  holes  in  p-type  silicon.  The  figures  also  show 
that  for  the  case  of  the  shallower  ionization  energies,  resistivity 
depends  more  strongly  on  temperature  for  the  lightly  doped  case  where 
lattice  scattering  is  dominant  and  become  less  temperature  dependent 
as  the  dopant  density  increases. 
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Theoretical  calculations  of  resistivity  vs  temperature  for 
indiuni-doped  silicon  with  dopant  density  as  a  parameter. 
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Figure  4.7.  Theoretical  calculations  of  resistivity  vs  dopant  density 
for  boron-doped  silicon  with  temperature  as  a  parameter. 


61 


10 


15 


1016 

_3 

Indium  Density  (cm  ) 


10 


17 


Figure  4.9.  Theoretical  calculations  of  resistivity  vs  dopant  density 
for  indium-doped  silicon  with  temperature  as  a  parameter. 


CHAPTER  V 


THE  HALL  FACTOR  IN  p-TYPE  SILICON 


5.1  Introduction 

The  most  direct  determination  of  the  mobility  is  by  the  Haynes- 

Schokley  drift  method,  wherein  the  drift  of  charge  carriers  in  a  known 

electric  field  is  measured.  However,  the  assumption,  made  when  these 

experiments  were  initiated,  that  the  drift  mobility  of  holes  as  minority 

carriers  in  an  n-type  sample  is  the  same  as  when  they  constitute  the 

majority  carriers,  is  invalid  in  view  of  carrier-carrier  scattering 

[59].  Also  the  experiment  can  succeed  only  if  the  lifetime  of  the 

minority  carriers  is  larger  than  the  transit  time.  For  this  reason, 

usually  Hall  mobilities  are  measured  instead.  The  Hall  mobility  is  the 

product  of  the  measured  conductivity  and  the  measured  Hall  coefficient. 

In  general  the  Hall  mobility  differs  from  the  conductivity  mobility  by 

a  factor  called  the  Hall  factor.  Determination  of  the  Hall  factor  may 

be  avoided  by  making  use  of  the  high  field  limit.  For  sufficiently 

high  magnetic  fields  several  simplifications  occur  in  the  magnetic  field 

dependence  of  the  Hall  coefficient.  In  the  high-field  limit  (when  the 
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product  of  mobility  and  magnetic  induction  becomes  greater  than  10  cm 
gauss/volt-sec  [60])  the  Hall  coefficient  is  simply  related  to  the 
carrier  concentration  by  [61] 
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and  thus  the  conductivity  mobility  and  the  Hall  mobility  are  equal. 
Although  the  high  field  limit  simplifies  use  of  the  Hall  mobility  con¬ 
siderably,  excessively  high  magnetic  fields  can  cause  problems  due  to 
the  quantization  of  the  hole  orbits  in  a  magnetic  field  [1].  The 
quantization  of  the  particle  motion  in  a  magnetic  field  will  create 
Landau  levels  within  the  band.  The  Landau  levels  will  modify  the  den¬ 
sity  of  states  in  the  valence  band  which  could  affect  the  interpreta¬ 
tion  of  experimental  data  [56].  Another  high  magnetic  field  effect  of 
importance  is  the  "magnetic  freeze  out"  which  occurs  with  the  stronger 
localization  of  bound  state  wave-functions  in  a  strong  magnetic  field 
[62].  Due  to  the  more  localized  charge  distribution,  the  Coulomb 
binding  energy  of  the  impurity  state  is  increased  so  that  at  a  fixed 
temperature  the  concentration  of  thermally  excited  charge  carriers  will 
be  smaller  and  the  Hall  coefficient  will  be  effectively  increased. 

Thus,  in  order  to  avoid  these  high  field  region  complications  and 
obtain  an  experimental  determination  of  the  value  of  conductivity 
mobility  in  the  low  field  limit,  it  is  necessary  to  have  an  accurate 
knowledge  of  the  Hall  factor  with  which  to  modify  measured  Hall 
mobilities.  Hall  measurements  are  routinely  used  to  experimentally 
determine  the  density  of  ionized  impurities  in  a  semiconductor  sample. 
This  determination  is  possible  only  if  an  accurate  value  of  the  Hall 
factor  for  the  particular  temperature  and  dopant  density  considered  is 
available. 

5.2  The  Hall  Factor 

The  Hall  and  conductivity  mobilities  are  related  by  the  Hall  factor 


as  follows: 


rH  -  yc 
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(5.2) 


For  nondegenerate,  spherically  symmetric  bands,  it  can  be  shown  that 
[59]  rH  a  1 ,  and  that 


rH  = 


<x2> 


<T> 


(5.3) 


In  general  most  previous  work  [5,59]  has  assumed  that  equation 

(5.3)  is  valid  in  the  case  of  p-type  silicon  and  thus  rH  will  vary 

between  1.18  (t(E)  E  ^  for  lattice  scattering)  and  1.93  (x(E)  <* 
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E  for  ionized  impurity  scattering),  if  hole-hole  scattering  is 
neglected,  and  will  approach  unity  for  the  degenerate  case.  This  theory 
does  not  allow  for  values  of  r^  less  than  one.  Experimental  evidence 
indicating  values  of  less  than  one  has  been  attributed  to  poor 
quality  of  the  measured  samples  [5].  Debye  and  Kohane  [63]  found  that 
the  measured  drift  mobility  for  holes  is  considerably  larger  than  the 
measured  Hall  mobility.  Values  of  r^  less  than  unity  were  also  reported 
by  Wolfstirn  [15]  for  the  case  of  gallium-doped  silicon.  More  recent 
experiments  [64]  show  that  a  value  of  r^  less  than  unity  is  necessary 
to  reconcile  differences  between  the  hole  concentration  measured  via 
Hall  coefficient  methods  and  that  inferred  from  dopant  densities 
determined  from  C-V  and  junction  breakdown  measurements.  The  usual 
assumption  made  is  to  let  r^  be  equal  to  one  and  thus  consider  the  Hall 
mobility  equal  to  the  conductivity  mobility.  Neglecting  the  Hall 
scattering  factor  alters  both  the  magnitude  and  temperature  dependence 
of  the  carrier  concentration  from  that  given  by  the  charge  balance 
equation.  In  fitting  data  to  the  charge  balance  equation,  both  thermal 
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carrier  concentration  and  dopant  impurity  activation  energy  are  over¬ 
estimated  by  the  assumption  of  unity  Hall  factor.  A  more  complete 
theoretical  treatment  of  the  Hall  factor  can  be  undertaken  by  consider¬ 
ing  the  nonparabolic  and  anisotropic  nature  of  the  valence  band  of 
silicon. 

Chapter  II  described  the  constant  energy  surfaces  as  warped  spheres. 
Warping  of  the  energy  surfaces  has  a  significant  effect  on  the  ratio  of 
Hall  to  conductivity  mobility.  When  the  bands  are  warped,  the  Hall 
factor  depends  on  the  degree  of  warping  as  well  as  the  scattering 
mechanism  [30]. 

The  Hall  mobility  is  the  product  of  the  ohmic  conductivity  and  the 
Hall  coefficient 

PH  =  °CRH  (5-4) 

In  the  low  field  limit  the  Hall  coefficient  for  a  nonparabolic, 
anisotropic  band  i  is  given  by  [37] 


R 


Hi 


(5.5) 


Thus  by  substituting  equations  (2.9)  and  (2.10)  into  equation  (5.5)  the 
Hall  coefficient  can  be  expressed  as 


R 


Hi 


?ie 


where 


(5.6) 


t 
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rHi  = 


fmCil 

(<Ti  "1 

(5.7) 


is  the  Hall  factor.  We  see  that  allowing  for  a  difference  between  the 
values  of  conductivity  and  Hall  effective  masses  due  to  the  anisotropic, 
nonparabolic  nature  of  the  band,  enables  us  to  separate  the  Hall  factor 
into  two  components:  the  mass  anisotropy  factor  given  by 


r 


Ai 


(5.8) 


and  the  scattering  factor  given  by 


rSi  = 


<Ti2> 


<Ti> 


(5.9) 


These  components  of  the  Hall  factor  will  be  considered  in  detail  in  the 
next  two  sections. 

5.3  The  Mass  Anisotropy  Factor 

Lax  and  Mavroides  [20]  have  derived  expressions  for  r^  based  on  the 
Dresselhaus  et  al  [28]  model  of  the  valence  band  of  germanium  and  sili¬ 
con.  Their  formulation  for  r^  acknowledges  the  anisotropy,  but  neglects 
the  nonparabolicity  of  the  bands.  In  general  it  is  found  that  r^  is  less 
than  unity  unless  the  scattering  anisotropy  becomes  extreme  [30].  In 
order  to  determine  the  variation  of  the  mass  anisotropy  factor  with 
changes  in  temperature  and  dopant  density  for  the  combined  valence  band 
of  silicon,  equation  (5.8)  was  evaluated  using  the  values  of  combined 
valence  band  effective  mass  obtained  from  equations  (2.28)  and  (2.29). 
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The  results  of  this  calculation  are  presented  in  Figures  5.1  and  5.2. 

These  figures  show  the  significant  contribution  of  the  mass  anisotropy 

factor  to  the  Hall  factor.  Since  the  influence  of  nonparabolicity  is 

reduced  in  degenerate  material  [25],  it  follows  as  shown  in  Figures  5.1 

and  5.2,  that  the  variation  of  r^  with  temperature  is  much  stronger  at 

low  dopant  densities,  since  it  is  in  this  dopant  density  range  that  the 

variation  of  effective  mass  with  temperature  is  the  strongest.  We  note 

that  the  mass  anisotropy  factor  is  less  than  unity  for  all  temperatures 

considered  in  this  work  once  the  dopant  density  increases  past 
15  -3 

6x10  cm  .  At  300  K,  r^  is  less  than  unity  even  for  dopant  densities 
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as  low  as  10  cm  . 

5.4  The  Scattering  Factor 

The  scattering  factor,  r<.,  depicted  in  Figures  5.3  and  5.4  as  a 
function  of  temperature  and  dopant  density,  does  not  follow  the  tradi¬ 
tionally  expected  variation  between  3tt/8  =  1.18  and  31  5tt/51  2  =  1.93  as 
the  dominant  scattering  mechanism  changes  from  lattice  to  ionized 
impurity  scattering.  Putley  [65]  has  noted  that  hole-hole  scattering 
can  modify  r^.  He  estimates  that  for  ionized  impurity  scattering,  r<. 
can  be  reduced  from  31  5tt/51  2  to  a  value  close  to  unity.  At  low  dopant 
densities  where  the  dominant  scattering  mechanism  is  acoustic  phonon 
scattering,  rs  varies  between  1.08  for  T  =  100,  to  1.24  for  T  =  400  K. 
The  deviation  from  the  traditionally  expected  value  of  r<.  =  1.18  is  due 
to  the  contributions  of  optical  phonon  modes  at  the  higher  temperatures. 
Hole-hole  collisions  also  affect  the  impurity  and  optical  phonon  scat¬ 
tering  contributions  so  they  become  significant  even  at  low  temperatures 
and  dopant  densities.  At  higher  values  of  dopant  density,  the  effects 
of  hole-hole  scattering  on  the  ionized  impurity  scattering  mechanism 
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Figure  5.1.  The  mass  anisotropy  factor  as  a  function  of  temperature 
for  various  impurity  dopant  densities. 
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Figure  5.4.  The  scattering  factor  r^  as  a  function  of  dopant  density  for  boron-doped 
silicon  with  temperature  as  a  parameter. 
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become  very  noticeable.  At  =  10  cm  ,  the  highest  value  of  r^  is 
1.29  for  T  =  100  K,  where  the  dopant  impurities  are  only  about  30  per¬ 
cent  ionized  [17].  At  higher  temperatures  where  the  percentage  of 
ionized  impurity  atoms  is  over  80  percent,  the  effects  of  hole-hole 
scattering  bring  r<.  from  its  traditionally  expected  value  of  1.93  to 
1.05  for  T  =  400  K. 

5.5  Hall  Mobility  and  Hall  Factor  in  the  Combined  Valence  Band 

Expressions  for  Hall  coefficient,  applicable  in  the  case  where 
holes  in  more  than  one  band  take  place  in  conduction,  are  given  by 
Putley  [66].  For  the  case  of  p-type  silicon,  assuming  no  compensation 
and  operation  in  the  low  field  region,  the  Hall  coefficient  is  given 
by  [66] 


R  =jjQ 

KH 


l  Rl^o, 


Hi  Ci 


3 

l  a, 

i=l 


Ci 


(5.10) 


By  substituting  equations  (2.9)  and  (5.6)  through  (5.10)  into  equation 
(5.4),  the  Hall  mobility  in  the  combined  valence  band  of  silicon  can 
be  expressed  by 
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The  conductivity  mobility  for  the  combined  valence  band  can  be  expressed 
by 


(5.12) 


Then  using  equations  (5.11)  and  (5.12)  we  can  express  the  Hall  factor  in 
terms  of  the  scattering  relocation  times  and  effective  masses  of  the 
individual  bands  by 


(5.13) 


Figures  5.5  and  5.6  summarize  the  results  of  equation  (5.13)  as  a 

function  of  temperature  and  dopant  density.  These  figures  show  that  the 

Hall  factor  ranges  theoretically  between  1.73  and  0.77  for  temperatures 

between  100  and  400  K  and  dopant  densities  between  10^4  and  1018  cm-8. 

For  temperatures  above  200  K,  r^  becomes  less  than  unity  for  dopant 
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densities  greater  than  5.5x10  cm  .  Figures  5.7  and  5.8  show  the 
theoretically  predicted  Hall  mobility  as  functions  of  temperature  and 
dopant  density.  These  two  figures  show  the  results  of  evaluating  equa¬ 
tion  (5.11)  with  the  aid  of  numerical  integration,  and  adjusting  the 
lattice  scattering  mobility  to  give  the  best  fit  to  values  of  conduc¬ 
tivity  mobility  deduced  from  resistivity  measurements. 
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Figure  5.5.  Theoretical  Hall  factor  vs  temperature  for  boron-doped 
silicon  with  dopant  density  as  a  parameter. 
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Figure  5.6.  Theoretical  Hall  factor  vs  dopant  density  for  boron-doped  silicon  with  temperature 
as  a  parameter. 
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Theoretical  Hall  mobility  as  a  function  of  temperature  for 
boron-doped  silicon  with  dopant  density  as  a  parameter. 


CHAPTER  VI 


EXPERIMENTAL  PROCEDURES 


6.1  Introduction 

Experimental  measurements  of  resistivity.  Hall  coefficient  and 

dopant  density  were  made  on  six  silicon  wafers,  four  doped  with  gallium, 

and  two  doped  with  indium,  in  dopant  densities  ranging  from  4.25xl015  to 
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3.46x10  cm  .  These  wafers  were  cut  from  crystals  grown  along  the 
<11 1 >  and  <100>  direction.  Additional  data  were  obtained  from  boron- 
doped  silicon  wafers  to  further  verify  the  adequacy  of  the  theory.  The 
data  were  obtained  from  test  patterns  NBS-4  [67]  fabricated  on  the 
silicon  wafers.  This  test  pattern  was  designed  at  the  National  Bureau 
of  Standards  primarily  for  use  in  the  evaluation  of  the  resistivity 
versus  dopant  density  relation  in  silicon.  Resistivity  measurements 
were  made  on  four-probe  square  array  resistors  and  collector  Hall  effect 
resistors,  while  the  net  dopant  density  in  the  specimens  was  determined 
by  the  junction  C-V  method  on  a  gated  base-collector  diode.  Mean  values 
of  resistivity,  dopant  density  and  Hall  coefficient  were  determined  by 
measuring  five  to  eight  selected  test  cells  with  a  standard  deviation 
in  resistivity  at  300  K  under  five  percent.  The  following  sections 

describe  the  test  sample  preparation  and  fabrication  procedure,  and  the 

% 

measurement  procedures. 
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6.2  Fabrication  Procedure 

The  overall  pattern -is  fabricated  on  a  square  silicon  chip  20 0  mils 
on  a  side  where  six  mask  levels  are  used  [68].  The  masks  were  used  in 
the  following  sequence:  base,  emitter,  base  contact,  gate  oxide,  contact 
and  metal.  Appropriate  cleaning  procedures  (see  Appendix  A)  precede  the 
diffusion  of  impurities,  and  a  negative  photoresist  process  was  used  in 
the  masking  steps.  The  base  mask  delineates  regions  whose  conductivity 
type  is  opposite  from  that  of  the  collector  substrate,  and  the  emitter 
mask  delineates  regions  whose  conductivity  type  is  the  same  as  that  of 
the  collector  substrate.  A  base  region  approximately  two  pm  deep  is 
diffused  into  the  background  material;  then  the  emitter  region  is  dif¬ 
fused  into  the  base  to  a  depth  of  approximately  one  pm.  The  base  con¬ 
tact  mask  is  used  to  open  windows  onto  the  base  region,  where  an  n+ 
diffusion  is  made  to  improve  ohmic  contact  to  the  base.  The  gate  oxide 
mask  delineates  regions  where  an  oxide  layer  of  closely  controlled 
thickness  is  grown  to  serve  as  a  gate  for  MOS  devices.  After  front-side 
metallization,  a  portion  of  the  wafer  was  separated.  This  section  was 
scribed  to  provide  the  Hall  effect  devices.  The  remainder  of  the  wafer 
was  then  metallized  on  the  backside  and  alloyed.  After  scribing,  the 
devices  were  mounted  on  T0-5  headers,  metal  contact  bonding  was  made, 
and  the  devices  were  encapsulated.  A  layer  of  ceramic  insulating  mate¬ 
rial  was  used  to  isolate  the  devices  from  contact  with  the  header. 
Resistivity  measurements  were  then  made  to  select  devices  for  use  in 
this  study. 
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6.3  Experimental  Measurements 

The  structures  used  to  evaluate  the  resistivity  of  the  bulk  mate¬ 
rial  are  the  Hall  effect  resistor  and  the  collector  four-probe  resistor 
[68].  The  four-probe  resistor  has  four  point  contacts  arranged  in  a 
square  array.  The  structure  (see  Appendix  A)  is  fabricated  by  diffusing 
a  base  over  a  large  area  except  at  the  four  point  contacts  which  are 
protected  from  the  base  diffusion  by  oxide  islands.  Emitters  are  dif¬ 
fused  at  these  points  in  order  to  make  low  resistance  contacts  to  the 
collector  material.  The  purpose  of  the  base  diffusion  is  to  eliminate 
surface  currents.  The  bulk  resistivity  is  determined  by  forcing  a 
current,  I,  between  two  adjacent  probes  and  measuring  the  voltage,  V, 
between  the  other  two  probes.  The  resistivity  of  the  material  is  deter¬ 
mined  from  [69] 
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where  S  is  the  probe  spacing  and  C  is  a  correction  factor  dependent  on 
the  ratio  of  probe  spacing  to  the  thickness  of  the  chip  [70].  This 
correction  factor  is  given  by 
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where  w  is  the  thickness  of  the  chip. 

The  collector  Hall  effect  resistor  is  a  four-terminal  resistor 
formed  in  a  square  chip  100  mils  on  a  side.  Contacts  are  formed  on  the 
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four  corners  by  an  emitter  diffusion  (see  Appendix  A).  The  resistivity 
is  calculated  from  [71] 
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where  w  is  the  thickness  of  the  chip,  and  V  is  the  voltage  difference 
between  nearest  neighbor  contacts  for  a  current,  I,  passed  between  the 
remaining  two  contacts.  The  TO-5  header  was  mounted  in  the  sample 
holder  of  an  Air  Products  and  Chemicals  AC-3L  CRYO-TIP  liquid  nitrogen 
system.  This  enabled  variation  of  the  sample  temperature  between  100 
and  350  K.  The  temperature  was  measured  by  a  chromel  vs  gold  with  0.07 
atomic  percent  iron  thermocouple.  For  temperatures  above  350  K,  the 
sample  was  placed  in  a  Stratham  Temperature  Test  Chamber. 

The  structure  used  for  the  Hall  coefficient  measurements  is  the 
collector  Hall  effect  resistor.  The  Hall  coefficient  is  calculated 
from 
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where  is  the  voltage  difference  measured  between  opposite  contacts 
for  a  current,  I,  passed  between  the  remaining  two  contacts,  and  B  is 
the  magnetic  field  density  perpendicular  to  the  plane  of  the  chip;  thus 
the  samples  are  oriented  so  that  the  magnetic  field  is  in  the  crystal 
growth  direction,  <11 1 >  for  the  gallium-doped  samples,  and  <100>  for 
the  indium-doped  samples.  The  Hall  mobility  is  determined  from 


02 


where  p  is  determined  from  resistivity  measurements  on  the  Hall  and 
four-point  structures.  The  magnetic  field  for  the  Hall  measurements 
was  provided  by  a  Varian  Associates  (V3703)  six-inch  electromagnet  with 
a  current  regulated  power  supply  (V-FR2503).  The  magnetic  field 
strength  was  monitored  by  a  Bell  620  gaussmeter  with  an  STB4-0402  probe 
with  a  stated  accuracy  of  0.1  percent.  Data  was  taken  over  a  tempera¬ 
ture  range  from  100  to  350  K.  The  current  used  in  the  resistivity  and 
Hall  coefficient  measurements  was  provided  by  a  Keithley  225  current 
source  capable  of  accuracy  within  ±0.5  percent  of  the  three-digit 
readout.  The  current  was  monitored  by  voltage  readings  across  precision 
resistors  connected  in  series  with  the  current  source.  These  resistors 
were  part  of  a  Dana-651  current  shunt  set,  accurate  to  within  0.01  per¬ 
cent.  Voltages  were  measured  with  a  Hewlett  Packard  3465A  digital 
multimeter  with  a  stated  accuracy  within  0.03  percent  of  the  readout. 
Resistivity  and  Hall  coefficient  measurements  were  made  in  accordance 
with  ASTM  standard  procedures  [72]. 

The  impurity  dopant  density  was  obtained  by  use  of  two  different 
structures:  an  M0S  capacitor,  and  a  base-collector  diode.  The  M0S 
capacitor  over  collector  consists  of  a  main  gate  which  is  surrounded 
by  a  field  plate  that  overlaps  a  channel  stop  which  also  serves  as  top 
side  collector  contact  [68].  This  structure  (see  Appendix  A)  is  used  to 
measure  the  collector  dopant  density  (N^  +  Nq)  from  the  high  frequency 
C-V  deep  depletion  method  [73].  The  collector  dopant  density  is  deter¬ 
mined  by  obtaining  a  dopant  profile  from  C-V  measurements  by  means  of 
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where  AV  is  an  incremental  change  in  the  gate  voltage,  and  the  measured 
capacitance  is  due  to  both  the  oxide  and  the  semiconductor.  A  self- 
consistent  check  was  made  on  the  measurements  of  collector  dopant  density 
by  using  the  base-collector  diode.  This  structure  (see  Appendix  A) 
consists  of  a  base  diffused  into  a  collector  and  a  metal  field  plate  to 
control  the  periphery.  The  field  plate  overlaps  both  the  base  and  a 
diffused  emitter  channel  stop  which  also  serves  as  topside  collector 
[68].  To  o.btain  a  correct  density  profile  the  field  plate  is  biased  at 
the  flat-band  potential  [74].  Capacitance-voltage  measurements  were 
taken  with  a  Princeton  Applied  Research  410  C-V  Plotter  and  a  Hewlett 
Packard  701 OA  X-Y  Recorder. 

From  each  silicon  wafer,  eight  four-probe  resistors,  eight  Hall 
resistors,  and  eight  capacitor-diode  chips  were  selected  for  encapsu¬ 
lation.  These  were  chosen  on  the  basis  of  low  leakage  currents  and  good 
contacts  at  the  metal  bonding  pads.  Measurements  were  made  on  each  of 
the  devices  and  data  from  the  five  to  eight  devices  closest  to  the  mean 
value  of  the  measurements  were  then  averaged.  In  this  manner  we  arrived 
at  representative  values  of  resistivity.  Hall  coefficient,  and  dopant 
density  for  each  sample.  The  results  of  these  measurements  and  compari¬ 
sons  with  the  theory  of  Chapters  III  through  V  are  presented  in  the 
next  chapter. 


CHAPTER  VII 


COMPARISON  OF  THEORETICAL  AND  EXPERIMENTAL  RESULTS 
7.1  Conduct! vi ty  Mobi 1 i ty 

Conductivity  mobility  was  evaluated  by  substituting  the  measured 

resistivities  into  equation  (4.7).  Figure  7.1  shows  the  hole  mobility 

plotted  as  a  function  of  hole  density  at  300  K  for  boron-doped  silicon. 

Curve  1  represents  the  theoretical  results  of  equation  (3.23),  and 

curve  2  was  reproduced  from  the  work  of  Wagner  [8].  Our  calculated 

values  are  within  six  percent  of  the  values  reported  by  Wagner  for 
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N ^  <  3x10  cm  .  For  higher  values  of  hole  density  our  calculated 
values  are  substantially  higher  than  those  of  Wagner.  As  previously 
explained  by  Li  [17],  this  discrepancy  is  due  to  Wagner's  assumption 
of  complete  ionization  of  boron  impurities.  This  assumption  is  valid 
only  at  low  dopant  densities  or  at  high  temperatures  where  full  ioniza¬ 
tion  of  boron  atoms  prevails.  The  theoretical  calculation  is  in 
excellent  agreement  (within  three  percent)  with  experimental  data 
reported  by  Thurber  et  al  [12].  This  gives  support  to  the  validity  of 
ionization  calculations  based  on  equations  (4.2)  through  (4.6). 

Mobility  values  reported  by  Horn  [13]  are  also  in  reasonable  agreement 
with  our  theoretical  results.  The  data  points  shown  in  Figure  7.1  were 
corrected  for  deionization  effects  via  equations  (4.2)  through  (4.6). 
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Figure  7.1.  Hole  mobility  vs  hole  density  for  boron-doped  silicon  at 
300  K.  Curve  1  is  the  theoretical  calculation  from 
equation  (3.23),  and  curve  2  is  reproduced  from  Wagner  [8]. 
Shown  also  are  the  experimental  results  of  this  work,  Horn 
[13],  and  Thurber  et  al  [12].  All  data  points  have  been 
corrected  for  deionization  effects. 
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7.2  Resistivity 

The  resistivity  vs  dopant  density  relationship  for  boron-doped 
silicon  at  300  K  is  shown  in  Figure  7.2.  The  solid  line  represents 
theoretical  calculations  using  equation  (4.7).  Wagner's  [8]  resistivity 
curve  and  the  theoretical  line  coincide  over  most  of  the  boron  density 
range.  Our  theoretical  calculations  agree  with  Wagner's  resistivity 
data  within  six  percent  over  the  entire  range  of  boron  densities  con¬ 
sidered  at  T  =  300  K.  Excellent  agreement  exists  between  our  experi¬ 
mental  data  and  the  theoretical  calculations  at  300  K.  Figure  7.2  also 
shows  excellent  agreement  between  our  theoretical  calculations  and  the 
data  of  Thurber  et  al  [12].  Good  agreement  was  obtained  with  the  data 
of  Thurber  and  Carpenter  [75]  where  total  boron  density  was  obtained  by 
the  nuclear  track  technique. 

Figure  7.3  shows  the  resistivity  of  gallium-  and  indium-doped  sili¬ 
con  as  a  function  of  total  dopant  density  for  T  =  300  K.  As  expected, 
because  of  the  deeper  ionization  energy  of  indium  as  compared  to  gallium, 
values  of  resistivity  for  gallium  doped  silicon  are  lower  than  values 
of  resistivity  for  indium-doped  silicon  at  the  same  total  dopant  density. 
Figure  7.3  does  not  show  this  at  low  dopant  densities  because  of  the 
assumed  values  of  background  boron  impurity  densities.  Excellent  agree¬ 
ment  was  obtained  between  our  experimental  data  and  that  obtained  from 
Wolfstirn  [15]  for  gallium-doped  silicon,  and  our  theoretical  calcula¬ 
tions  at  T  =  300  K.  Data  obtained  from  the  two  indium-doped  samples 
showed  good  agreement  with  the  theoretical  calculations,  but  the  same 
was  not  true  for  the  data  of  Schroder  et  al  [64],  and  Backenstoss  [16]. 
As  seen  in  Figure  7.3,  for  each  value  of  measured  resistivity,  Schroder 
et  al  [64]  report  two  different  values  of  measured  indium  density.  The 
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Figure  7.2.  Resistivity  vs  dopant  density  for  boron-doped  silicon  at 
300  K.  Shown  are  the  theoretical  calculations  from 
equation  (4.7)  and  the  experimental  results  of  this  work, 
Thurber  and  Carpenter  [75],  and  Thurber  et  al  [12]. 
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lower  value  of  indium  density  was  obtained  by  C-V  and  junction  breakdown 
methods,  while  the  higher  value  was  obtained  by  Hall  measurements  and 
curve  fitting.  Because  of  uncertainties  in  the  value  of  the  Hall  scat¬ 
tering  factor,  Schroder  et  al  consider  the  lower  value  of  density  more 
reliable.  Note  that  our  theoretical  calculation  falls  between  the  two 
values  of  dopant  density  reported  by  Schroder  et  al  [64].  Values  of 
resistivity  of  indium-doped  silicon  reported  by  Backenstoss  [16]  are 
about  25  percent  higher  than  our  calculated  values.  The  work  of 
Backenstoss  [16],  however,  was  done  in  the  high  doping  region  where 
dopant  densities  approach  the  limit  of  solid  solubility.  Backenstoss 
found  that  for  dopant  densities  greater  than  4x10^  cm”^  there  was  a 
considerable  amount  of  indium  precipitation.  Thus  it  is  possible  that 
part  of  the  discrepancy  between  our  theoretical  calculations  and  the  data 
of  Backenstoss  is  due  to  the  low  solid  solubility  limit  of  indium  in 
silicon.  Recent  theoretical  results  of  Sclar  [14]  for  In-doped  silicon 
also  agree  very  closely  with  our  theoretical  calculations  at  300  K. 

To  find  out  the  adequacy  of  our  theoretical  model  for  temperatures 
other  than  300  K,  we  compared  the  calculated  and  measured  values  of 
resistivity  for  silicon  samples  doped  with  boron,  gallium  and  indium 
for  temperatures  ranging  from  100  to  400  K.  Figure  7.4  shows  the  com¬ 
parison  between  the  theoretical  and  measured  resistivities  for  boron- 
doped  silicon.  Except  for  a  couple  of  data  points,  agreement  between 
the  theoretical  and  measured  values  was  within  8  percent  over  the  entire 
range  of  temperatures.  Figures  7.5  and  7.6  show  the  comparison  between 
theoretical  and  measured  resistivities  for  gallium-  and  indium-doped 
silicon  respectively.  Agreement  here  was  not  as  good  as  in  the  boron- 
doped  case,  but  except  for  a  couple  of  data  points,  agreement  between 
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Figure  7.4.  Resistivity  vs  temperature  for  the  boron-doped  silicon 

samples.  Solid  lines  are  the  theoretical  calculations  and 
solid  dots  are  the  experimental  data. 
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Figure  7.6.  Resistivity  vs  temperature  for  the  indium-doped  silicon 
samples. 


theory  and  experiment  was  within  10  percent.  For  the  indium-doped 
samples,  the  largest  discrepancies  occurred  at  low  temperatures.  This 
is  suggestive  of  some  degree  of  compensation  in  the  samples.  An  experi¬ 
mental  estimate  of  percentage  of  compensation  was  not  made  for  any  of 
the  silicon  samples  studied  in  this  work. 

7.3  Hall  Mobility 

Figures  7.7  through  7.12  show  the  experimentally  determined  Hall 
mobility  for  silicon  slices  doped  with  gallium  and  indium.  Agreement 
between  theory  and  experiment  is  good  for  temperatures  above  150  K  for 
the  gallium-doped  samples.  However,  for  temperatures  of  150  K  and  below 
some  points  differ  from  the  theoretical  prediction  by  as  much  as  38 
percent.  The  measured  Hall  mobility  of  the  indium-doped  samples  is  well 
within  15  percent  of  the  calculated  value  except  for  one  data  point  at 
T  =  100  K.  Our  control  over  the  temperature  of  the  samples  was  better 
for  the  T  >  200  K  range,  but  the  main  source  of  the  discrepancy  between 
theoretical  and  experimental  values  of  Hall  mobility  is  the  Hall  factor. 
This  will  be  discussed  in  the  next  section.  The  magnetic  field  was 
rated  accurate  to  within  1  percent.  Data  points  representing  a  Hall 
mobility-dopant  density  pair  are  estimated  to  have  a  total  error  of 
around  8  percent. 

7.4  Hall  Factor 

The  Hall  factor  in  the  case  of  p-type  silicon,  is  plotted  as  a 
function  of  dopant  density  for  T  =  300  K  in  Figure  7.13.  With  the  ex¬ 
ception  of  the  points  deduced  from  Morin  and  Maita's  [5]  data  for 
boron-doped  silicon,  agreement  between  the  Hall  factor-dopant  density 
data  and  the  theoretical  calculation  is  within  15  percent  for 
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Figure  7.8.  Hall  mobility  vs  temperature  for  gallium-doped  sample. 
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Figure  7.9.  Hall  mobility  vs  temperature  for  gallium-doped  sampl 
Nfl  =  1.26xl017  cm’3. 


Figure  7.10.  Hall  mobility  vs  temperature  for  gallium-doped  sample. 
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Figure  7.11.  Hall  mobility  vs  temperature  for  indium-doped  sample. 
Na  =  4. 64x1 Q1 6  cm"3. 
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Figure  7.12.  Hall  mobility  vs  temperature  for  indium-doped  sample. 
Na  =  6. 44x1 01 6  cm'3. 


Hall  factor  vs  dopant  density  for  p-type  silicon  at  300  K.  Shown  are  the 
theoretical  calculations  from  equation  (5.13)  and  the  experimental  results 
of  this  work,  Wolfstirn  [15],  Morin  and  Maita  [5],  and  Schroder  et  al  [64] 
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>  5x10  cm  .  While  in  general,  agreement  between  theory  and 
experiment  in  this  region  of  dopant  density  is  only  within  15  percent, 
the  bulk  of  the  experimental  data  is  within  10  percent  of  the  theoreti¬ 
cal  prediction.  The  samples  used  by  Morin  and  Maita  [5]  showed  a  con¬ 
siderable  degree  of  compensation;  in  their  lower-doped  samples,  the 
compensation  was  as  high  as  20  percent.  This  fact  may  have  contributed 
to  the  low  values  of  measured  Hall  mobility,  and  the  large  discrepancy 
between  these  data  points  and  the  calculation.  This  model  neglects  the 
effects  of  compensation,  and  the  combined  presence  of  acceptor  impurities 
of  varying  ionization  energies.  The  overall  effect  of  adding  impurities 
of  both  signs  is  one  of  increased  ionization  of  the  excess  or  ND  [15]. 
It  would  be  necessary  to  know  the  percentages  of  compensation  of  the 
experimentally  measured  crystals  to  accurately  determine  the  adequacy 
of  the  theory  at  low  dopant  densities.  An  experimental  estimate  of 
percentage  of  compensation  was  not  made  for  any  of  the  silicon  samples 
studied  in  this  work.  Long  [38]  has  noted  that  the  low  measurements  of 
Hall  mobility  for  p-type  silicon  may  not  be  due  entirely  to  compensation 

and  the  quality  of  the  crystals.  For  reasonably  pure  silicon  samples 
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(p  =  35  ohm-cm,  =  4.4x10  cm"  )  Long  [30]  obtained  Hall  mobilities 

2 

between  360  and  390  cm  /volt-sec.  Hall  mobilities  as  great  as  450 
2 

cm  /volt-sec  have  also  been  reported  [76].  While  higher  than  the  Hall 
mobilities  of  Morin  and  Maita,  these  measurements  [38,76]  still  indicate 
a  value  of  Hall  factor  for  low-doped  p-type  silicon  at  300  K  less  than 
unity.  It  is  still  doubtful  that  a  Hall  mobility  smaller  than  a  con¬ 
ductivity  mobility  at  300  K  is  really  an  intrinsic  property  of  p-type 
silicon  [38],  It  is  possible  that  the  Hall  factor  may  be  greater  than 
unity  in  a  crystal  of  exceptionally  high  perfection.  However,  calcu- 


102 


lations  involving  an  exact  model  of  the  valence  band  of  silicon  [77] 
indicate  that  the  anisotropy  of  the  bands,  while  not  important  in  the 
interpretation  of  mobility  and  resistivity  data,  plays  an  important 
role  in  the  evaluation  of  the  Hall  factor.  Our  model  averages  out  the 
anisotropy  of  the  bands  and  thus  our  values  of  the  mass  anisotropy 
factor,  r^,  are  larger  than  an  exact  model  would  predict.  A  calculation 
based  on  the  exact  model,  applicable  over  the  entire  range  of  interest 
of  temperatures  and  dopant  densities,  is  beyond  the  scope  of  this  work. 
At  the  higher  dopant  densities  there  is  fair  agreement  between  the 
theoretical  prediction  and  values  of  Hall  factor  deduced  from  our 
experimental  work  and  the  data  of  Wolfstirn  [15].  Data  points  from 
Wolfstirn's  work  showing  ionization  energies  greater  than  0.065  eV  were 
not  included  in  Figure  7.13  because  of  the  high  degree  of  compensation 
of  these  samples. 


CHAPTER  VIII 


SUMMARY  AND  CONCLUSIONS 

Theoretical  expressions  have  been  derived  to  compute  the  hole 
mobility,  resistivity  and  Hall  factor  as  functions  of  dopant  density 
and  temperature  for  silicon  doped  with  boron,  gallium  and  indium.  The 
valence  band  of  silicon  was  represented  by  a  three-band  model  which 
takes  into  account  the  nonparabolic  nature  of  the  bands.  This  attri¬ 
bute  of  the  valence  band  is  included  in  the  effective  mass  calculations. 
Contributions  from  scattering  by  acoustical  and  optical  phonons, 
ionized  impurities  and  neutral  impurities  were  considered  in  the  calcu¬ 
lation  of  average  relaxation  time.  In  addition,  our  model  also  takes 
into  account  the  effect  of  hole-hole  scattering  on  both  lattice  and 
ionized  impurity  scattering  relaxation  times,  and  the  effect  of  inter¬ 
band  transitions  on  the  acoustic  phonon  relaxation  time.  Thus  the  model 
developed  in  this  study  represents  a  more  complete  theoretical  descrip¬ 
tion  of  the  conductivity  mobility.  Hall  mobility.  Hall  factor  and 
resistivity  than  previous  theoretical  models  have  acknowledged. 

As  stated  in  Chapter  II,  the  model  used  for  the  valence  band  of 
silicon  provides  values  of  density-of-states  effective  mass  which  lead 
to  values  of  intrinsic  carrier  density  which  are  well  within  the  limits 
of  experimental  error.  The  greatest  deficiency  this  model  has  over  an 
exact  calculation  is  that  the  anisotropy  of  the  bands  has  been  averaged 
out.  Thus  while  the  temperature  dependence  of  effective  mass  derived 
from  the  exact  and  approximate  models  of  the  valence  band  is  similar, 
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the  magnitudes  of  the  masses  may  be  substantially  different.  This  dif¬ 
ference  is  not  relevant  in  the  study  of  mobility  and  resistivity  in 
p-type  silicon,  but  becomes  important  in  the  determination  of  the  Hall 
factor.  Since  direct  experimental  verification  of  values  of  conductivity 
and  Hall  effective  masses  is  not  possible,  the  only  way  to  assess  the 
value  of  effective  mass  calculations  is  by  using  the  theory  in  the 
development  of  directly  measurable  properties  such  as  resistivity  and 
Hall  mobility. 

The  resistivity  analysis  for  the  boron-,  gallium-,  and  indium-doped 
silicon  samples  showed  agreement  between  experimental  and  theoretical 
results  within  10  percent  over  the  entire  range  of  temperature, 

100  s  T  <  400  K.  Note  that  best  agreement  between  theory  and  experiment 
was  obtained  for  boron-doped  samples,  followed  by  gallium-  and  indium- 
doped  samples.  This  may  be  due  to  the  fact  that  we  neglect  the  compen¬ 
sation  effect  and  the  possible  dependence  of  ionization  energy  on  dopant 
density  in  the  theoretical  calculations  for  gallium-  and  indium-doped 
samples.  An  experimental  estimate  of  degree  of  compensation  was  not 
made  for  any  of  the  silicon  samples  studied  here.  Data  points  repre¬ 
senting  a  resistivity-dopant  density  pair  are  estimated  to  have  a  total 
error  of  around  6  percent. 

A  comparison  between  our  calculated  mobility  values  with  those  of 

Wagner's  [8]  data  on  boron-doped  silicon  at  300  K  shows  that  agreement 
.  17  -3 

is  within  ±6  percent  for  <  3x10  cm  .  Discrepancies  at  higher 
dopant  densities  can  be  eliminated  if  the  effect  of  deionization  of 
boron  impurities  were  included  in  Wagner's  calculations  [8].  Excellent 
agreement  was  found  between  our  theoretical  calculations  of  mobility  in 
boron-doped  silicon  and  the  data  of  Thurber  et  al  [12]  at  300  K.  We 
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have  obtained  excellent  agreement  between  theoretical  values  of  resis¬ 
tivity  and  our  experimental  data  for  boron-,  gallium-,  and  indium-doped 
silicon  at  300  K;  our  theoretical  calculations  also  agreed  with  the 
resistivity  data  by  Thurber  et  al  [12]  for  boron-doped  silicon,  and 
Wolfstirn  [15]  for  gallium-doped  silicon. 

As  shown  in  Chapter  V,  however,  this  theoretical  model  does  not 
have  the  same  kind  of  success  in  describing  the  Hall  mobility  and  Hall 

factor  for  p-type  silicon  for  low  dopant  densities,  and  low  temperatures. 

1 6  -3 

For  dopant  densities  lower  than  5x10  cm"  ,  this  model  predicts  a  value 
of  Hall  factor  much  greater  than  has  been  experimentally  determined 
[5,15,38],  A  discrepancy  like  this  for  the  case  of  p-type  germanium 
was  eliminated  as  the  quality  of  germanium  crystals  improved.  Thus  it 
has  been  assumed  that  low  values  of  Hall  mobility  for  low-doped  p-type 
germanium  at  T  =  300  K  were  caused  by  compensation.  However,  compensa¬ 
tion  alone  may  not  account  for  low  values  of  Hall  mobility  in  the  case 
of  p-type  silicon  [38].  Recent  studies  by  Nakagawa  and  Zukotynski  [77] 
indicate  that  the  use  of  the  exact  model  in  the  development  of  the  Hall 
factor  formulation  yields  results  which  agree  in  general  with  experimen¬ 
tal  data  for  the  case  of  p-type  silicon.  Experimental  data  for  the 
case  of  p-type  germanium  does  not  agree  with  the  theoretical  results 
of  Nakagawa  and  Zukotynski  [77]. 

From  this  study,  we  have  found  that  the  theoretical  expressions 
derived  in  this  work  are  adequate  for  mobility  and  resistivity  calcula¬ 
tions  for  p-type  silicon  in  the  temperature  range  100  <  T  <  400  K,  and 
the  dopant  density  range  10^  <  <  10^  cm  The  theoretical  formu¬ 

lation  is  also  adequate  for  description  of  the  Hall  factor  and  Hall 

16  _  *3 

mobility  for  dopant  densities  above  5x10  cm  and  temperatures  above 
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200  K.  The  failure  of  this  model  in  predicting  Hall  factor  and  Hall 
mobility  at  low  temperatures  and  dopant  densities  is  due  to  the  omission 
of  the  proper  anisotropy  formulation  in  the  model  of  the  valence  band 
of  silicon.  This  weakness  in  the  model  would  be  remedied  by  use  of  the 
exact  model  in  a  manner  similar  to  that  of  Nakagawa  and  Zukotynski 
[11,77]. 

Further  improvements  to  this  model  would  include  the  exact  formu¬ 
lation  of  the  valence  band  of  silicon  on  Kane's  [27]  model,  as  prescribed 
by  Nakagawa  and  Zukotynski  [11,77].  Thus  the  proper  nonparabolicity 
and  anisotropy  would  be  included.  For  better  comparison  with  experi¬ 
mental  results,  the  variation  of  ionization  energy  with  dopant  density 
should  be  included  in  the  calculation  of  hole  density  for  the  deeper 
impurities.  Also  the  effects  of  compensation  and  the  presence  of  other 
p-type  impurities  in  silicon  samples  must  be  considered.  It  would  also 

[  be  of  great  benefit  to  extend  this  study  into  the  heavy  doped  region. 
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APPENDIX  A 


FABRICATION  PROCEDURES  AND  TEST  STRUCTURES 

This  appendix  contains  a  list  of  the  fabrication  procedures 
followed  to  generate  the  test  structures  measured  in  this  study.  These 
test  structures  were  part  of  NBS-4  [67]  test  pattern,  and  include  a 
planar  four-probe  collector  resistor,  a  100  mil  square  Hall  effect 
device,  a  gated  base-collector  diode,  and  a  gated  MOS  capacitor  over 
collector  structure. 

A.l  Initial  Si 02  Masking: 

A. la  Initial  Clean  Up: 

(1)  Ultrasonic  clean  in  hot  DI  water  with  small  amount  of 
100  Tritonex  solution  for  10  minutes. 

(2)  Rinse  in  running  DI  water  for  5  minutes. 

(3)  Place  in  solution  of  1 NH^OH: 1 H2O2 : 2H2O  for  20  minutes 
at  50  degrees  C. 

(4)  Rinse  in  DI  water  for  5  minutes. 

(5)  Dip  in  10  percent  HF  for  10  seconds. 

(6)  Rinse  in  DI  water  for  5  minutes. 

(7)  Place  in  solution  of  1 HC1 : 1 HgOg : 2H2O  for  20  minutes  at 
50  degrees  C. 

(8)  Rinse  in  DI  water  for  10  minutes  and  spin  dry  in  Ng. 

A. lb  Initial  Oxidation  -  350  nm  at  1100  degrees  C 

(1)  Push-in  5  minutes  Ng  at  1000  cc/minute 
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(2) 

Dry  02 

5  minutes 

02  at  800  cc/minute 

(3) 

Wet  02 

35  minutes 

02  at  800  cc/minute 

(4) 

Dry  02 

5  minutes 

02  at  800  cc/minute 

(5) 

n2 

15  minutes 

N2  at  1000  cc/minute 

(6) 

Pull-out 

5  minutes 

N2  at  1000  cc/minute 

A. 2  Phosphorus  Base  Diffusion: 

A. 2a  Photoresist  (PR)  Application  with  Base  Mask  (NBS-4-1AB): 

(1)  Bake  at  200  degrees  C  for  30  minutes  to  completely  dry 
the  wafer. 

(2)  Apply  Waycoat  200  negative  PR  (do  this  on  both  sides  of 
the  wafer,  apply  and  spin  PR  on  back  side  first). 

(3)  Spin  at  5000  RPM  for  20  seconds. 

(4)  Prebake  in  65-degree  C  oven  for  20  minutes. 

(5)  Align  wafer  and  mask  and  expose  for  4  seconds  under  UV 
light.  (Both  sides  of  the  wafer  must  be  exposed.) 

(6)  Develop  for  20  seconds  in  spray  of  Waycoat  developer 
(undiluted),  then  clean  off  developer  with  15-second 
spray  of  Butyl  Acetate. 

(7)  Dry  with  and  inspect  under  the  microscope. 

(8)  Post  bake  at  130  degrees  C  for  25  minutes  in  or  air 
circulating  oven. 

A. 2b  Base  Window  Etch: 

(1)  Etch  in  buffered  HF  for  2  minutes,  rinse  in  DI  water, 
dry  and  inspect.  Etch  for  15  seconds  more  to  see  if 
the  pattern  changes  color. 
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(2)  When  certain  that  all  the  oxide  has  been  removed,  place 
wafer  in' hot  (90  degrees  C)  J-100  PR  remover  for  5 
minutes. 

(3)  Quench/rinse  in  methanol  for  5  minutes,  then  rinse  with 
methanol  from  squirt  bottle. 

(4)  Rinse  in  DI  water  for  5  minutes. 

(5)  Clean  up  -  same  as  above  A. la  except  for  ultrasonic 
clean. 

A. 2c  Phosphorus  Base  Diffusion  at  875  degrees  C: 

(1)  Push-in  5  minutes  ^  at  1250  cc/minute 

(2)  Phosphorus  16  minutes  N2  at  1250  cc/minute 

02  at  100  cc/minute 
N2  through  source 
bubbler  at  10  cc/ 


minute 

T(source)  =  4  to  7 
degrees  C 

(3)  Flush-out  2  minutes  O2  at  800  cc/minute 

Ng  at  1250  cc/minute 

(4)  Pull-out  5  minutes  N2  at  1250  cc/minute 

A. 2d  Removal  of  Phosphosilicate  Glass: 

(1)  Dip  in  10  percent  HF  for  5  seconds. 

(2)  Rinse  In  DI  water  for  5  minutes  and  spin  dry  in  N?. 

A.2e  Drive-in  Diffusion  and  Base  Oxide  -  350  nm  at  1100  degrees  C, 

base  depth  approximately  2pm: 

(1)  Push-in  5  minutes  N2  at  1000  cc/minute 

(2)  Dry  02  10  minutes  02  at  800  cc/minute 


45  minutes 


02  at  800  cc/minute 


(3)  Wet  09 

(4)  Dry  0^  10  minutes  0^  at  800  cc/minute 

(5)  Drive-in  70  minutes  at  1000  cc/minute 

A. 3  Boron  Emitter  Diffusion: 

A. 3a  PR  with  Emitter  Mask  NBS-4-2AB  -  Same  as  A. 2a  (front  side  PR 
only),  with  the  following  change:  Before  applying  Waycoat 
PR,  coat  surface  of  the  wafer  with  silazane  (5  parts 
Hexamethyldisilizane:  95  parts  Xylene),  and  spin  at  4000 
RPM  for  20  seconds. 

A. 3b  Emitter  Window  Etch: 

(1)  Etch  in  buffered  HF  for  2.5  minutes,  rinse  in  DI  water, 
dry  and  inspect.  Etch  for  15  seconds  more  to  see  if 
the  pattern  changes  color. 

(2)  When  certain  that  all  oxide  has  been  removed  (oxide 
islands  must  be  gone),  strip  PR  and  clean  up  same  as 
A. 2b  (5). 

A. 3c  Boron  Nitride  Diffusion  at  1050  degrees  C: 

(1)  Load  wafers  facing  source. 

(2)  Push-in  5  minutes  N2  at  500  cc/minute 

(3)  Pre-dep  30  minutes  N2  at  500  cc/minute 

02  at  25  cc/minute 

(4)  Oxide  10  minutes  02  at  1000  cc/minute 

(5)  Pull-out  5  minutes  N2  at  500  cc/minute 

A. 3d  Removal  of  Borosilicate  Glass: 

(1)  Dip  in  10  percent  HF  for  10  seconds. 

(2)  Rinse  in  DI  water  for  5  minutes,  spin  dry  in  N2> 


A.3e  Boron  Drive-in  and  Emitter  Oxide,  T  =  925  degrees  C: 


(1) 

Push-in 

5 

ini  nutes 

N2 

at 

1000  cc/minute 

(2) 

Dry  02 

5 

minutes 

°2 

at 

1000  cc/minute 

(3) 

Wet  02 

65 

minutes 

°2 

at 

1000  cc/minute 

(4) 

Dry  02 

5 

minutes 

°2 

at 

1000  cc/minute 

(5) 

Passivation 

5 

minutes 

n2 

at 

1000  cc/minute 

(6) 

Pul  1 -out 

5 

minutes 

A. 4  Base  Contact  Diffusion: 

A. 4a  PR  with  Contact  Diffusion  Mask  NBS-4-5AB  -  Same  as  A. 2a  (PR 
on  both  sides  -  back  side  first).  Back  side  must  be  exposed. 
A. 4b  Base  Contact  Window  Etcn: 

(1)  Etch  in  buffered  HF  for  2  minutes,  45  seconds.  Rinse 
dry,  check  for  oxide  removal. 

(2)  When  all  oxide  has  been  removed  strip  PR  and  clean  up 
same  as  A. 2b  (5). 

A. 4c  Phosphorus  Base  Contact  Diffusion  at  875  degrees  C: 

(1 )  Push-in  5  minutes 

(2)  Pre-dep  10  minutes 


(3)  Flush-out  2  minutes 


at  1250  cc/minute 
at  1250  cc/minute 
0^  at  100  cc/minute 
through  source 
bubbler  at  10  cc/ 
minute 

T( source)  =  4  to  7 
degrees  C 

O2  at  800  cc/minute 
at  1250  cc/minute 
N~  at  1250  cc/minute 


(4)  Pull-out 


5  minutes 
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A.4d  Removal  of  Phosphosi 1 icate  Glass: 

(1)  Dip  in  10  percent  HF  for  5  seconds. 

(2)  Rinse  in  01  water  for  5  minutes  and  spin  dry  in  N^. 
A.4e  Reoxidation  at  95  degrees  C: 


(1) 

Push-in 

5  minutes 

N2  at  1000  cc/minute 

(2) 

Dry  0^ 

10  minutes 

02  at  1000  cc/minute 

(3) 

Wet  02 

25  minutes 

02  at  1000  cc/minute 

(4) 

Ory  02 

5  minutes 

02  at  1000  cc/minute 

(5) 

Passivation 

5  minutes 

N2  at  1000  cc/minute 

(6) 

Pull-out 

5  minutes 

MOS  Gate  Oxide: 

A. 5a 

PR  with  Gate  Oxide  Mask  NBS-4-6AB 

-  Same  as  A. 2a  (PR  on  front 

side  only,  use  silizane). 

A. 5b 

Etch 

in  Buffered  HF  for  3  minutes. 

rinse,  dry  and  check  for 

oxide  removal. 

A. 5c 

When  all  oxide  has 

been  removed  strip  PR  and  clean  up  same  as 

A. 2b  (5). 

A.5d 

Gate  Oxide  Growth  ■ 

o 

■  12GQ0A  at  950  degrees  C: 

(1) 

Push-in 

5  minutes 

N2  at  1000  cc/minute 

(2) 

Dry  02 

10  minutes 

02  at  800  cc/minute 

(3) 

Wet  02 

40  minutes 

02  at  800  cc/minute 

(4) 

Dry  02 

10  minutes 

02  at  800  cc/minute 

(5) 

Passivation 

10  minutes 

N2  at  1000  cc/minute 

(6) 

Pull-out 

5  minutes 

A. 6  Contact  Windows: 

A. 6a  PR  with  Contact  Mask  NBS-4-3AB.  Same  as  A. 2a  (front  side 
only,  use  silizane). 
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A. 6b  Contact  Window  Etch: 

(1)  Etch  in  Suffered  HF  for  3  minutes,  check  for  complete 
removal  of  oxide. 

(2)  Strip  PR  and  clean  up  same  as  A. 2b. 

A. 7  Front  Side  Metallization: 

A. 7a  In  the  alloy  furnace  set  at  450  degrees  C,  5  minutes  in  N2; 

15  minutes  in  forming  gas,  5  minutes  in  N2>  Gas  flow  =  800 
cc/minute  in  all  cases. 

A. 7b  Evaporate  800  nm  aluminum  over  the  front  surface  of  the  wafer. 
A. 7c  PR  with  Metal  Mask  NBS-4-4AB.  Same  as  A. 2a  (front  side  only, 
use  silizane). 

A.7d  Etch  in  A1  etchant  of  20  H^PO^H^IHN  heated  to  50  degrees 
C  until  the  etching  is  complete  (about  2  minutes). 

A.7e  Rinse  in  01  water  for  10  minutes. 

A.7f  Strip  PR  in  J-100  and  methyl  as  prescribed  in  A. 2b  (3  minutes 

in  80  degrees  C  J-100). 

A. 8  Back  Side  Contact: 

A. 8a  Dry  for  30  minutes  in  200-degree  C  oven. 

A. 8b  Spin  PR  on  front  side  as  prescribed  in  A. 2a  (no  mask). 

Expose  to  U.V.  light  for  4  seconds. 

A. 8c  Hard  bake  at  130  degrees  C  for  30  minutes. 

A.8d  Dip  in  10  percent  HF  until  oxide  has  been  removed  from  the 
backside  of  the  wafer. 

A.8e  Remove  PR  in  J-100  and  methyl  as  prescribed  in  A. 2b  (do  not 
use  cleaning  procedure). 

A.8f  Scribe  the  wafer  to  separate  the  Hall  devices. 
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A.8g  Evaporate  400  nm  aluminum  over  the  back  side  of  chips. 
A. 9  Post  Evaporation  Micro  Alloy.  Load  chips  in  alloy  furnace  at 
500  degrees  C  for  15  minutes.  N2  set  for  a  flow  rate  of 
approximately  100  cc/minute. 


Figure  A. 3.  Gated  base-collector  diode.  Diameter  =  17  mils  (431.8  ym) 
NBS  -  4.8  [67], 


APPENDIX  B 


EXPERIMENTAL  SETUP  AND  DATA 

This  appendix  consists  of  a  block  diagram  which  shows  the  experi 
mental  setup  and  equipment  used,  and  six  tables  which  show  the  data 
taken  on  the  test  structures.  In  addition  to  measured  values  of  resi 
tivity,  dopant  density  and  Hall  coefficient,.  Tables  B-l  through  B-6 
also  contain  values  of  calculated  Hall  mobility. 
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Keithley-225 


Dana-651 


#589876 


Figure  B.l.  Experimental  Set-up.  For  Hall  effect  measurements  the  mag¬ 
netic  field  was  generated  by  a  Varian  Associates  6  in. 
Electromagnet  System  #544772.  Capacitance  measurements 
were  taken  with  a  Princeton  Applied  Research-410  C-V  Plot¬ 
ter  #601922  and  a  HP-7010A  X-Y  Recorder  #602659. 
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Table  B-l.  Sample  SGA  16-2  (Ga-Doped  Si,  =  4. 25x1 01 5  cm"3) 


Temperature 

T(K) 

Resistivity 

(ft-cm) 

Hall  Coefficient 

2 

R^(cm  /coulomb) 

Hall  Mobility 

2 

PH(cm  /v-sec) 

100 

1 .700 

9381 

5518 

120 

1.185 

140 

1 .066 

150 

1 .063 

2239 

2106 

160 

1.097 

200 

1.434 

1443 

1006 

250 

2.230 

1433 

643 

300 

3.300 

1263 

383 

350 

4.410 

1244 

282 

400 

5.580 

t 


122 


Table  B-2.  Sample  SGA  16-3  (Ga-Doped  Si,  =  4.09xl016  cm'3) 


Temperature 

T(K) 

Resistivity 

(ft-cm) 

Hall  Coefficient 

R^(cm  /coulomb) 

Hall  Mobility 

o 

uH(cm  /v-sec) 

100 

0.660 

2325 

3523 

150 

0.308 

413 

1341 

160 

0.295 

180 

0.286 

190 

0.288 

200 

0.294 

206 

701 

250 

0.361 

163 

452 

300 

0.469 

149 

318 

350 

0.593 

147 

248 

400 


0.735 
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Table  B-3.  Sample  SGA  17-2  (Ga-Doped  Si,  =  1 .26x1 01 7  cm-3) 


Temperature 

T(K) 

Resistivity 

(ft- cm) 

Hall  Coefficient 

R^Ccm  /coulomb) 

Hall  Mobility 

o 

yH(cm  /v-sec) 

100 

0.464 

1019 

2196 

150 

0.204 

199 

975 

180 

0.176 

200 

0.172 

88 

512 

220 

0.175 

250 

0.188 

65.7 

349 

300 

0.216 

57.8 

268 

350 

0.253 

52.1 

206 

400 


0.294 


124 


Table  B-4.  Sample  SGA  18-2  (Ga-Doped  Si,  NA  =  3.46xl017  cm"3) 


Temperature 

T(K) 

Resistivity 

(il-cm) 

Hall  Coefficient 

3 

R^(cm  /coulomb) 

Hall  Mobility 

2 

uH(cm  /v- sec) 

100 

0.550 

605 

1100 

150 

0.163 

101 

620 

200 

0.116 

46.2 

398 

230 

0.111 

240 

o.m 

250 

0.112 

30.7 

274 

300 

0.121 

25.2 

208 

350 

0.135 

21.8- 

161 

400 


0.152 
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Table  B-5.  Sample  SIN  16-2  (In-Doped  Si,  =  4.64x10^  cm 


Temperature 

T(K) 

Resistivity 

(ft-cm) 

Hall  Coefficient 

R^(cm  /coulomb) 

Hall  Mobility 

o 

PH(cm  /v-sec) 

100 

21.150 

71571 

3384 

150 

4.73 

9899 

2093 

200 

1 .840 

1855 

1008 

250 

1.131 

602 

532 

300 

0.916 

307 

335 

350 

0.877 

207 

236 

400 

0.913 

!  26 


Table  B-6.  Sample  SIN  17-1  (In-Doped  Si,  NA  =  6.44xl016  cm’3) 


Temperature 

T(K) 

Resistivity 

(ft-cm) 

Hall  Coefficient 

3 

R^(cm  /coulomb) 

Hall  Mobility 
uH(cm2/v-sec) 

100 

25.95 

58958 

2272 

150 

4.65 

7268 

1563 

200 

1.564 

1345 

860 

250 

0.952 

458 

481 

300 

0.766 

237 

309 

350 

0.721 

158 

219 

400 

0.740 

f 


APPENDIX  C 


COMPUTER  PROGRAM 

This  appendix  contains  the  FORTRAN  computer  program  used  to 
generate  values  of  conductivity  and  Hall  effective  mass,  conductivity 
and  Hall  mobility,  and  Hall  factor  as  functions  of  temperature  and 
acceptor  impurity  density  in  silicon.  Besides  the  main  program,  two 
subroutines  are  also  listed.  The  first  of  these,  DENSTY,  calculates 
the  hole  density  by  iterating  the  Fermi  energy  level  in  the  charge 
balance  equation  as  explained  in  Chapter  IV;  the  second,  INTBND, 
computes  the  coefficients  (functions  L..  and  T . .  of  equations  (3.7) 

I  J  I  J 

through  (3.9))  used  in  the  calculations  of  acoustic  phonon  scattering 
based  on  the  work  of  Bir  et  al  [43]. 
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IOC 

20C  CONDUCTIVITY  EFFECTIVE  MASS.  HALL  EFFECTIVE  MASS. 

30C  CONDUCTIVITY  MOBILITY.  HALL  M03ILITY.AND 

40C  HALL  FACTOR  AS  FUNCTIONS  OF  TEMPERATURE  ANO 

45C  ACCEPTOR  IMPURITY  DENSITY  IN  SILICON 

SOC  DOPED  WITH  BORON 

60C 

70  REAL  MDI .MD2.MD3. MD.NN.NAI .NA.LIMI 

60  REAL  MI.M2.M3.MIT.M2T.MJT.K.K0.M0.KT 

90  REAL  MH.MH1 .MH2. MH3 .MC .MCI . MC2.MC3 

100  REAL  MOB  1 .M0B2.M0B3 .MOB .MOBH 

110  REAL  MOBZ 

12  OC 

130  DIMENSION  C R DP( 1 01  .PC TG < 50 )  . T ( 1 0 » 

140  DIMENSION  NN { 50 ) • N A  I { 50 ) . PP( SO ) * N A ( 50) 

150  DIMENSION  MD 1  (1  0  )  •  MD2  (  l  0  )  .  MD3  ( l  0  )  •  MO  (1  0  ) 

160C 

170  F5 (X)=(X-C3)**1 .5*EXP(C3-X 1 

160  F6(X)=X**1.5*EXP(-X) 

190  F7(X)=(X-C2)**1.5*EX»(-X) 

20  OC 

210  DATA  K, A, 0 «C /8. 62E-5 .4. 27.0.63.4.93/ 

220  OATA  DELTA. MO. KO/O. 044, 9. 10956E-31. 1.38062E-23/ 

230  DATA  HBAR.H.PI/1.05459E-34.6.6252E-34.3.14159/ 

240  DATA  EO.ES/3. 854 185E- 12. 11. 7/ 

250  DATA  TAUO .MW/6.960E-1 0. 0.244/ 

260  DATA  0/1.602195-19/ 

270  2>*TA  THETA.  M  1. M2.  M3/735.0.  0.537,  0.  153.0.234/ 

280C 

29 OC  READ  IN  TEMPERATURES  AND 
300C  DENSITY  OF  STATES  EFFECTIVE  MASSES 
31 OC  READ  IN  THE  H-H  SCATTERING  CORRECTION 
31  SC  FACTOR  TO  OPTICAL  PHONON  SCATTERING 
320C 

330  RE  AD (5.10)<T(J)i  3=1.10  1 

340*10  FORMAT ( 10F6.0I 

350  READIS.l 1 1 (MDI (LI *L=1 .10 1 

360  READ(  5.11)(MD2(L)  »L=1  .101 

370  READ(5. 11  )(M03(L )*L=1*  101 

380  RE  AO (5.1 1 1 ( MD (L) »L=1 .101 

390*11  FORMAT! 10(F 7.5. 1X1 1 

400  RE AD (5 . 1 2 ) ( CROP (L 1 *L= I  .  10) 

410*12  FORMAT!  1  0(F6.4.1  X)  1 

420C 

430C  READ  IN  DOPANT  CONCENTRATIONS 
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440C 

430  REAO  (5.13) (NAIL). L=t, 40) 

460413  F  OR MAT (5(E8.2*1X) ) 

47  OC 

480C  COMPUTE  F (GAMMA ) • S  AND  ASSPCT  AT  FD  CONSTANTS 
49 OC  FOR  THE  DETERMINATION  OF  EFFECTIVE  MASSES 
SOOC 

310  8P  =  SORT(B**2*4(C**2./6.0n 

520  APB=A+9P 

525  AMB=A-BP 

530  GM=C**2./(2.*BP*AMB> 

540  GP=(C**2./( 2«*BP*AP6) )*(- 1. ) 

550  F0M=1 .04.05 *GM4*0 1 635 *GM** 2 . ♦ .0 0 0 9 08*GM**3 • 

560  F2M=1. 4. 0I667+SM+. 041 35 9*S M* *2 . 4 . 00090679* GM**3 • 4 

570  1*0  0091 959*GM  **4.04.00002106*GM**5. 

580  e2P— 1  *  4 • 0 1 667*GP4 «041359*GP**24.  000906  79*GP**3 ♦ 

590  1 *0 009 1959*GP **44. 000021 06*GP**5 

600  F3M*l .-0.0 1667* GM40. 01 7956* GM**2.- 

602  60*  0069857*GM* *3*  4*001 261  *GM**4 

610  F3 P= 1 • —0*0  I  66  7*GP  40  *01 7956*GP**2 

612  C-0.0069857*GP**34.0012510*GP**4 

620C 

630C  COMPUTE  SCATTERING  LIFETIME  CONSTANTS 
64  OC 

650  TOPI l=WW*THETA/TAUO 

660  T0P22=T0P11 

670  TOP33=TOPll 

680  T I A® 0*4 2* /SORT! MO) 

690  TIB=0**2./(4 .0*PI*ES*E0 )**2. 

700  TIC=TIA*TIB/( KO*  * 1 • 5*4 »  5015BE— 7) 

710  GA«7.539822E-S*MO/0**2. 

720  GB=« (ES*E0*4.0*PI )/H) *( K0**2./H> 

730  GC=GA*GB 

740  TNE= 1*1 36E- 1 9/ES  **2 • 

750  TNF®8.  88576 E6*HBAR** 2*0/ (M0**1 *5*SORT(KO I) 

760C 

77 OC  THIS  LOOP  CHANGES  THE  TEMPERATURE 
78  OC 

790  DO  200  J*  1*8 

800  KT=K*T(JJ 

810  C2=DELTA/(3.*KTI 

820  C3*DELTA/KT 

830  C4»THETA/TI J> 

840  C5*C34C4 

850  C6«C3— C4 
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660 

670 

660 

890 

900 

910 

920 

930*30 

931 

940 

950 

960 

970*40 
90OC 
99 0C 
1000C 
1010 
1020 
1030 
1040C 
1050C 
1060C 
1070C 
1060 
1090 
1100 
1101 
1110 
1120 
1130 
1140 
1150 
1160 
1170 
1160 
1190 
1200 
1210 
1220 
1230 
1240 
1250 
1260 
1270 
1280 


R8=(M01 ( J)/MD< J) )**1.5 
R9=CMD2( J)/MD( J)  >4*1.5 
RIO* (MD3<  J)/MD(  J  >  >**l  .5*CXP(-C3> 

EG*1. 205-2. 8E-44T( J) 

CR  OGAP*  (  1.205/204*0.66667 

M3T=M3*CBDGAP 

WRITE (6.30) 

FORMAT i 1H1////8X. ‘CONDUCT! VI  TY  AND  HALL  EFFECTIVE 
C-MASS.  CONDUCTIVITY  AND  HALL  MOBIL  ITY. ‘/8X,  ‘AND 
6HALL  FACTOR  AS  FUNCTIONS  0=  DOPANT  CONCENTRATION, 
t  BORON  DOPED  SILICON*///) 

WRITE(6.40)T(  J) 

FORMAT! /40X. «T*‘.P4. 0////1 

COMPUTE  INTERBAND  SCATTERING  PARAMETERS 

Y1*MD1 ( J> 

Y2*MD2l  J> 

CALL  TNTBNDIYl.  Y2.  C  t  1,  C22.C  12.  C21.CI  12> 

INTRODUCE  TEMPERATURE  AND  MASS  PARAMETERS 

TO  THE  CALCULATION  OF  SCATTERING  RELAXATION  TIMES 

AC11T=(C11/TAU0>*MD  KJI  **1 .5*TI  J)  **  1.  5 
AC22T=CC22/TAU0 )*MD2t J  ) ** I .5*7 C J  > 4*  1. 5 
AC 12T=(C 12/TAUO) 4M01 < J > 4 4 1 . 5*T <  U) *4 1. 54 
&IMD2CJ  >/MDl ( J) >**2.S 
AC21T*(C21/TAU0)4MD2( J >4*1 .5*T(J  >441.5 
T  AC33= 1 • 0/T  AUO 
AC33T*T AC334T (J)441.5 
P0IST*1./«SXPCC4)-1  .) 

PDIST1=P0IST*I. 

0P1 1  TasTOPI  1  *SQRT(T(  J>  >*MD1  f  J  >441 .5 
0P22T  »T0P22*S0RT(  T(  J)  >  *  MD  2(  Jl  *  4  l  .  5 
OP33T=TOP33*SORTI T( J  >>  4M3T4*  1.5 
Till  T=TI  C/ (T (J)4*l .5*  MD1 ( J  >440.5) 

TI22T  =TIC/iT<  J ) 44 1 , 5*  MD2( J>**0.5) 

T I33T =T IC/IT { J ) 4*  1  • 54M3T  4*0.5) 

B 1T=GC*  T(J)**2. 4MD1 ( J> 

B2T=GC*TCJ  )  442.4MD2I  J) 

B3T=GC*TO)  442.4M3T 
TNG— TNE *MD  (  J)/(  K04TC  J>  > 

TNHI *TNF/ (5QPT(T< J) )*MDI(  J  >441.5) 
TNH2=TNF/(S0RT(TC  J) ) 4  MD2 ( J  >  *41  .5  > 

TNH3=TNF/( SORT! T( J> )*MD3IJ >441. S> 


v 

i. . 


v 
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1290C 

1300C  COMMUTE  SCATTERING  LIFETIMES 

131 OC  THIS  LOOP  ALSO  CHANGES  THE  CONCENTRATION 

1320C 

1330  00  100  Lsl.40.2 

1340C 

13S0C  CALL  SUBROUTINE  TO  CALCULATE  IONIZED  OOPANT  DENSITIES 
1360  C 

1370  QI=TCJ1 

13S0  0  2=NA(  L  I 

1390  Q3*M0(J) 

1400  CALL  DENSTYIQl  .02.03.04,05.06.07,08  ) 

1410  N A I { L )— Q  4 

1420  PP( L ) =05 

1430  NN( L) =Q6 

1440  PCTG<L»  =  Q7 

1450  ETA=08 

1460C 

1470C  COMPUTE  CORRECTION  FACTORS  FOR  HOLE-HOLE  SCATTERING 
1460C 

1490  IF(NA(L)  *LT.  1.0EI5)GO  TO  44 

1500  IF(NAIL)  .LT.  3.0E17IGO  TO  46 

1510  C  RAC  =  0«  86 

1520  GO  TO  48 

1530*44  CRAC— 1 • 0 

1540  GO  TO  48 

1550*46  CRAC=1.0004-{4.0I3378F-19*NA(_»> 

1560*48  CRI=(NAI (L)/PP(L) ) • C 1 . -EX P(-PP ( L )/NA I (L ) )) 

1570C 

1580C  INTROOUCE  CONCENTRAT rONS  AND  H-H  SCATTERING 
1581C  CORRECTIONS  INTO  RELAXATION  TIME  CALCULATIONS 
1600C 

1610  A Cl  1 S=ACl 1 T/CRAC 

1620  AC22S=AC22T/CRAC 

1630  AC1 2S= AC  1 2T /CR A  C 

1640  AC21 S=AC21 T/CRAC 

£650  AC33S=AC33T/CRAC 

1660  OP l IS=0PI IT/ CROP  CO ) 

1670  0P22S  =0P22T  /CROPCJ) 

1680  OP33S=OP33T/CROP( Jl 

1690  G1*B1T/PPCL» 

1700  G2=B2T/PP<L) 

1710  G3*83T/PP(L  ) 

1720  GBI=AL0G(G1  M.l-Gl/IGITI.  > 

1730  G82sAL0G(G2+U  )-G2/(G24  1.  > 


r 
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1740 
1750 
1760 
1770 
1780 
1790 
1800 
18  10C 

1820C  COMPUTE  AVERAGE  RELAXATION  TIMES  BY  SIMP  SON’S  RULE 
1830C  INTEGRATION.  FOR  BAND  2.  100  ITERATIONS  TO  GET  TO 
1832C  THE  •KNEE*  OF  THE  BAND*  AND  200  TO  CONTINUE  BEYOND 
184 OC  THE  • KNEE* 

18S0C 

I860  LlMl  =  0.02/{K*TC  Jl) 

1870  E  NO— SO 

1880  N  1* 100 

1890  N2=200 

1900  N3=300 

1910  .  W=L  I  M  J  /(  2  *N1) 

1920  Z=(END-LIMI >/I2  *N2I 

1930  SUM| =0. 0 

1940  SUM  2=0.0 

1950  SUM3=0.0 

1960  SUM4-=  0.  O 

1970  SUM5=0.0 

1980  SUM6=0.0 

1990  SUM  8=0.0 

2000  SUM9=0 .0 

2010  SUM1 0=0.0 

2020  SUM  11-0.0 

2030  SUM12=0.0 

2040  SUM1 3=0.0 

2050  SUM  1 9=0 • 0  , 

2060  SUM20=0.0 

2070  SUM 2 3=0. 0 

2080  TLA33  =0.0 

2090  X3=0 . 0 

2100  X  4=  0. 0 

2110  X  5=  0 . 0 

2120  X6=0. 0 

2130  T 133=0. 0 

2140  V=0.0 

21S0C 

2160C  SET  UP  SIMPSON • S  RULE  INCREMENTS 
2170C 


GB3=AL0G ( G3  ♦  1  • 1 -G3/ (G3*l . I 
Till S=TI 11T4NAI  ( L ) *  GB l /CRI 
TI22S  =T I22T  AN  AI  <  L  >  *GB  2/C  R  I 
TI33S=TI33T*NAI CL) AG83/CR I 
TNI 1=TNH1*NNCL) 
TNI2=TNH2*NN1L) 

TNI  3*TNH3ANN4L) 


I 
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2180  00  70  KK=1 . N3 

2190  IF(KK  .GT.  HUGO  TO  101 

2200  GO  TO  102 

22108101  XL*V*2.*(  KK-N1  -1  )*Z 

2220  XR*V«-2.*(KK-N1  )*Z 

2230  XM«(XL«-XR>/2. 

2240  W»Z 

2250  GO  TO  103 

2260*102  XL»2.*(KK-1)4W 

2270  XR»2.*KK*W 

2280  XM*(XL»XR>/2. 

2290  IF(KK  .CO.  Nl)V=XR 

2300*103  IF(XL  «C0.  0. 0 ) XL*XM/ l 000. 0 

2310  00  65  JJ*1 .3 

2320  IF  (JJ  .CO.  1)  X»XL 

2330  IF  (JJ  .CO.  2)  XsXB 

2340  IF  (JJ  *E0.  3)  X=XM 

235  OC 

2360C  ACOUSTIC  PHONON  SCATTERING  RELAXATION  TIMES 
2370C 

2380  XI  *  SORT ( X  1 

2390  TLAll  — AC1 1 S*X  1 

2400  TLA  1 2  sAC125*Xt 

2410  TLA22  «AC22S*X1 

2420  TLA21  *AC2 1 S*X1 

2430  TA33  »AC33S*X1 

2440  I F( XL  .LT.  C31G0  TO  54 

2450  TLA33  =AC33S* SORT ( X  -C3 ) 

2460*54  CONTINUE 

2470C 

2480C  OPTICAL  PHONON  SCATTERING  RELAXATION  TIMES 
2490C 

2500  I F ( XL  -LT.  C4 ) GO  TO  57 

2510  X3  »SORT(X  -C4) 

2520*57  X4  *SORT(X  *C4) 

2530  TLOll  =0P1 1 S*(PDIST1 4X3  ♦PDISTFX4  ) 

2540  TL33  *OP33S* ( PO I  ST 1* X3  +P3IST4X4  ) 

2550  IF  (XL  .LT.  C5 ) GO  TO  58 

2560  X5  *  SORT ( X  -C5> 

2570*58  IF (XL  .LT.  C6IG0  TO  59 

2580  X6  =SQRT(X  -C6 ) 

2590*59  TL033  =0P33S* ( PO I ST1 * X5  fP0IST*X6  I 

2600  TL022*0P22S*  (PDIST  l*X3 *-P0 1 ST*X4  ) 

261  OC 

2620C  IONIZED  IMPURITY  SCATTERING  RELAXATION  TIMES 
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2630C 
2640 
2650 
2660 
267  0 
2660 
2690*63 
2700C 
27 10C 
2720C 
2730 
2740 
2750 
2760 
277  0C 
2760C 
2790C 
2800C 
2810 
2820 
2830 
2840 
28S0C 
2860 
2870 
2880 
2890 
2900 
2910 
2920 
293001 
2940 
2950 
2960 
2970 
2980 
299002 
3000 
3010 
3020 
3030 
3040 
305003 
3060 
3070 


X8  *1.0/X  4*1.5 
Till  =TI11S*X6 
TI33P  =TI33S*X8 
IF!  XL  .LT.  C3)G0  TO  63 
T 133  -T  I33S/IX  -C3>**1.5 
TI22  *T 122  S  *X8 

NEUTRAL  IMPURITY  SCATTERING  RELAXATION  TIME 

X9=X  l  +TNG/X  1 
TNl*TNIl/X9 
TN2=TNI2/X9 
T  N3*TN 13/ X 9 

T0T4L  RELAXATION  TIMES  IN  EACH  BAND 
INTRABAND  SCATTERING  TIMES 

T 1  sl.O/CTLAll  +TL011  4TI11  4TN1) 

T2  *1#0/(TLA22  *TL022  *7122  *T N2 ) 

T3  *  1 • 0/C  TA 33  4TL33  +TI33P  4TN3J 
TAU3  =  1 .0/(  TLA33  4TL033  *TI33  «-TN3) 

D  *  1 .0— T 1  *T 2  /ITLA12  *TH2l  » 

CR1  =1.0+(T2  4MD1 ( J ))/ (TLA12  44D2IJ)) 

CR2  =  1.0+CTt  4MD2I J)  )/( TLA 2 1  4  MD 1 ( J ) ) 

T  AU 1  *CR1  *T 1  /D 
TAU2  -CR2  *T2  /D 

TZ»1 .0/IR8/TAU14R9/TAU24R10/TAU3I 
GO  TO  (1 .2,3). JJ 
TAU1L=TAUI 
TAU2L=TAU2 
T3L»T3 
TAU3L=TAU3 
T2L=TZ 
GO  TO  4 
TAU1R»TAU1 
TAU2R=TA U2 
T3R*T3 
TAU3R=TAU3 
TZR»T 2 
GO  TO  4 
TAUIM=TAUI 
TAU2M=TAU2 
T3M*T3 


K 

i 

b 


f*r 


k 


-v  Strain 

■ 
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3080  TAU3M=TAU3 

3090  TZM=TZ 

31 0004  CONTINUE 

3! 10 #65  CONTINUE 

3120C 

3130C  TOTAL  LIFETIMES  IN  EACH  BAND  SQUARED 
3140C 

3150  T 1 SQL=T  AUI L**2 • C 

3160  T1SQM=TAU1M**2. 0 

3170  T 1SQR=TAU1R**2.0 

3180  T2SOL=TAU2L**2.0 

3190  T2S0M=TAU2M**2. 0 

3200  T2SQR=TAU2R**2.0 

3210  T33S0L=T3L**2. 0 

3220  T33SQM=T3N  **2.0 

3230  T33  SQR=T3R**2 .0 

3240  T3S0L  =  TAU3L**2. 0 

3250  T 3SQM=T  AU3M  **2.0 

3260  T3SQR=TAU3R**2.0 

3270  TS0ZL=TZL**2.0 

3280  TSQZM=T ZM*  *  2.0 

3290  TSQZR=TZR**2.0 

3300C 

331 OC  SUMMATIONS  FOR  THE  INTEGRALS 
332 OC 

3330  SUM  1* SUM  1  ♦(  W/3.  )  *C  T 3L*F6(  XL  )*4.  *T3M*>=6(  XM  )  * 

3335  6T3R*F6(XR)) 

3340  SUM 3= SUM 3*  C W /3. 1*{ TAU2L*F6( XL) * 

3345  64  .*T AU2M*F6  <XM )  *T AU2R*F6 (XR ) ) 

3350  SUM6=SUM6*(  W/3.  )*(  TAU1  L*F6( XL)  ♦ 

3355  64  .*T AUI M*F6 (X* ) *T AU IR*F6( XR ) ) 

3360  SUM8=SUM8* ( W/3 . ) * ( T33S  QL*F6 (XL ) * 

3365  64.*T33SOM*F6(XM)*T33SQR*F6(XR) ) 

3370  SUM10=SUM10*(W/3. )*(T2S0L*F6( XL )* 

3375  C4.*T2SQM*F6(  XM )  *T2  SQR* F6 ( XR )  ) 

3380  SUM 13= SUM  13* (W/3. )*( TISQL*F6( XL)* 

3385  t4.*TlSQM*F6(XM)*TlSQR*F6(XR) ) 

3390  SUM  1 9=  SUM 1 9* ( W/3. ) * ( TZ L*  =  6 I  XL) * 

3395  64.*TZM*F6(XM)*TZR*F6(XR)> 

3400  SUM20=SUM2  0* (W/3 . ) * ( T5 QZL*F6 (XL ) * 

3405  64.*TS0ZM*F6( XM)*TSQZR*F6(XR) ) 

3410  I F ( XL  .LT.  C3)G0  TO  66 

3420  SUM2=SUM2*< W/3. ) * ( TAU3 L* F5 ( XL) + 

3425  64 . *T  AU3M4F5 ( XM ) *TAU3R*F5( XR )  ) 

3430  SUM9=SUM9*(W/3. )* ( T3SQL *F5 ( XL ) ♦ 
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3435 

3440 

3450466 

346  0 

3470 

3475 

3480 

3485 

3490 

3500*67 

3510 

3515 

3520 

3525 

3530*68 

3540*70 

3550C 

3560C 

3S70C 

3580 

3590 

3600 

3610 

3620 

3630 

3640 

3650 

3660 

3670 

3680 

3690 

3700 

37  IOC 

3720C 

3730C 

374  OC 

3750 

3760 

3765 

3770 

3780 

3790 

3795 

3800 

3810 


&4.AT3S0M4F5{ XM) ♦TSSQRTFSI  XR>  ) 
SUM23=SUM23MX/3.  >*(  F5(  X_  >  4  4.4F51  XM  >«-F5l  XR>  > 
CONTINUE 

IFIKK  .GT •  N1JG0  TO  67 
SUM4sSUM44(W/3.  > * < T AU 2L *F6 ( XL  >♦ 

£4  •  4  TAU2M4F 6  (  XM  )  4-TAU2  R*  F6  (  XR  >  ) 

SUM1  Is  SUM  1  l+U/3.  )*(  T2SGL*F6( XL >♦ 
64»4T2SQM*F5(XM)  4-T2S0R4F6  C  X  R  )  ) 

GO  TO  68 
CONT INUE 

SUM5  =  SUM5M  W/3.  >  *  (  T  AU2  L*F7  (  XL  > 

64  .*TAU2M*F  7( XM ) +TAU2R4F  7< XR»  ) 

SUM1  2  =  SUM l  2  +  It / 3»  1*(T2S0l*F7(XL)«- 
64. 4T2SQM4F7I XM) ♦ T2 SQR* F7 { X R > ) 

CONTINUE 

CONTINUE 

AVERAGE  LIFETIMES 

SUM21sl. 32934 
TAVl=SUM6/SUM21 
TAV2  =SUM3/SUM21 
T AV 3sSUM2/SUM23 
TAV Z— SUM1 9/SUM2 1 
T  SO V 1=  SUM1 3/SUM21 
TSQV2  =  SUM 10/SUM21 
TS0V3=SUM9/SUM23 
T  SO  VZ- SUM2  0/SUM2 1 
TPAT  t=TSQV 1 /TAV 14*2.0 
TPAT2*TS0V2/TAV2**2.0 
TRAT3=TSQV3/TAV3*42.0 
TRATZ«TSOVZ/TAVZ**2 .0 

THE  INTEGRALS  HAVE  BEEN  CALCULATES 
NOW  COMPUTE  THE  MASS 

R 7s SUM 1 /SUM  7 

R3=SUM3/t ( F2P» SUM4/SORT CAPB» )♦ 

C(  F2M6SUM5/SQRTI  AMB  > I > 

RS=MOI I J)**l .5/<S0RTf AM6J4F3M ) 

R6*M03< J)4* 1 • 54 SUMS* EXP ( C3 ) / (SORT ( A  )*SUM9 ) 
R7*SUM 1 0/ ( F3P*S0RT< APB )*SUM1 1* 

•F3M*  SORT ( AMB  ) 4 SUM l 2 ) 

MC1*M01C  J)**l.5*SQRTCAMB) /F2M 
MC2*MD2(  J  >441.5  *R  3 
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3620  MC3=MD3!  J) * *1 • S*S0RT  C A >*CXP! C3  > *>R2 

3830  M  HI  =  5QR T  !  R  5  J 

3840  MH2®S0RT !MD2!J)**1 »5*R7) 

3850  MH3*  SORT !  R6) 

3860C 

387 OC  CALCULATE  MC.MH.HALL  FACTOR.  AND  MOBILITY 
3880C 

3890  R 1 1=T AV 1 / ! T AVZ *MC 1 ) 

3900  »12®TAV2/!  TAVZAMC2  > 

3910  R  13®TAV3/!TAVZ*MC3) 

3920  R14*TS0V1/!TS0VZ*MHI 4*2.0 ) 

3930  R 15*TS0V2/( TSO VZ*MH2**2.0» 

3940  R16*TSQV3/<TS0VZ4MH3**2.0) 

3950  MC  *!•/!  (Rll*R8)MR12*R9)MR134R10)  ) 

3960  MH  ®l  ./SORT<  (R84P14)  «■!  R9*R15)MR1  0*R16)  ) 

3970  M081 * 1 .0  E4*T AY 14Q/IMC1 *M0> 

3980  MOB 2® 1 • OF  4*  TAV2  *0/(MC2*M01 

3990  M0B3® 1.0E4*TAV34Q/CMC3*M0» 

4000  M0BZ=1.0E4*TAVZ*a/!MC*M0) 

4010  W0B®RQ*M0BH-P9*M092*R1  0*M093 

4020  RES® 1.0/ ( Q*MOB*NA 1 1  L  )  ) 

4030  RHZ®TRATZ*( MC/MHJ442.0 

4040  MOBH®MOB*RHZ 

4240C 

4250C  WRITE  OUT  COMPUTEO  VALUES  OF  MOBILITY.  MASS 
42S5C  RESISTIVITY.  ETC. 

4260C 

4270  WRITE  I  6. 85) NAIL) . NAI (L)  .PP ( L ) , NN( L ) .PCT GIL ) 

4280*85  FORM  AT  CSX.  *N  A®  •  .  E  1  0 . 4,  3X.  ^=*.*10.4, 

4285  £3  X.  *  PP=» ,El 0.4 ,3X . *NN=* .E10.4.3X.  • N AI /N A® • . F6. 47) 

4300  WRITE! 6.92) MOlf  J)  *  MD2 I J )  • MD3  C  J )  «  MD I J) 

4310*92  FORM  AT  (  5X  •  •  MD1  =  •  .  F7. 4,  5X,  •  M3  2=  •  .  F  7.  4.  5X. 

4315  £• MD3  =  » .F7.4.5X.* M0®» ,F7.4/) 

4330  WR I TE ( 6.94) MC 1 . MC  2. MC 3. MC 

4340*94  F0RMATI5X,*  MCI =• .F7 .4 . 5X. • MC2* • . F7.4. 5X. 

4350  t»MC3-=.F7.4,5X.»  MC  =  *  .F7.4/) 

4360  WRITEI6. 96 )MH1 .MH2. MH3 . MH 

4370*96  =  OR  MA  T I  5  X  » • MH1 =•  .F7.4.5X.* MH2  =  *  . F7.4.5X. 

4380  £*  MH3®  •  «F  7.  4  »  5X  .  •  MH®  *.c7.4/) 

4390  WRITE  16 .88  JTAV1.TAV2.TAV3.T  AVZ 

4400*88  e0RMAT!5X, »TAVl=» ,E 1 0. 4 ,5 X , • TA V2 * * .E10.4.5X. 

4410  £»TAV3® ».E10. 4. 5X. • TA V® • . E l 0. 4/) 

4420  WRI TE 16 .89 ) TS0V1 . TS0V2 . TSQV3 . TS QV Z 

4430*89  FORMAT! 5X.»TSQVI  =  *,E10 .4 .5 X. • TSQV2** • El  0. 4 *5X • 

4440  £»TSQV3=* .E10.4.5X, *TSQV® • , E 10. 4/) 


i 


138 


31  MAY  1979 


4450 

4460*93 

4470 

4510 

4S20486 

4530 

4540 

4550*98 

4560 

4590 

4600*99 

4610 

4620C 

463  0C 

4640*100 

4650*200 

4660 

4670*300 

4660 

4690 


WRI TF(6.93) TRAT1 .TRAT2 .TRAT3 .TRATZ 

FORMAT! SX.  *TRATl=*.F6.3.5X.*TRAT2  =  *.ir6.3.5X. 

C*  TRAT3=*  .F6.3.5X.*TRATZ=*.F6.3/) 

WRITEC6.86)C11 .C22.C12  .C21  „  . 

FORMAT ! 5X«  •  Cl  1=*.F6.3.5X«*C  22= • .F6«3.5X.*C12=*. 
1F6. 3.5X.*  C2  1=*  »F6«3/) 

WR I TE ( 6. 98) MP3  1 .MOB 2 . M^B3 .MOP Z 
FORMAT  (  5X  •  •  MOB  l  =  *  »  F8«  2. 5X  «  *'4032=*.F8.2.  5X. 
fc»  M0B3=*  .F9. 2  .5  X . • MOBZ=*  .  FB.2/> 

JR1TEI6. 99)MCB.MOBH.RHZ,RcS 
FORMAT !5X .* COND  MOB= • . F3 • 2 . 3X . • HALL  MOB= * • F 8 • 2. 3X. 
1 • HAUL  FACTOR**  ,F  7„  4  .3  X  ,  •  PS  S  I STI  VI  TV  =•  .F8.3////V  1 


CONTINUE 
CONTINUE 
WRI TE (6.300) 
FORMAT! IH1) 
STOP 
END 
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4790  SUBROUTINE  OENSTY  <01 . 02 .03  •  04  .  Q5  •  06.  Q7 .08  ) 

4700C 

4710C  SUBROUTINE  TO  CALCULATE  IONIZED  AMD  NEUTRAL 

4715C  IMPURITY  DENSITY 

4720C 

4730C  NV=EFFECTIVE  DENSITY  OF  BAND  STATES 

4740C  EA=DOPANT  IONIZATION  ENERGY 

4750C  NA=DCPANT  IMPURITY  CONCENTRATION 

476 OC  £G=ENSRGY  GAP 

4770C  EF=FERMI  ENERGY 

4780C 

4800  REAL  NV . NA . K ,KT. MO 

4810  REAL  NAI.NN 

4820  DATA  K/8.6173E-5/ 

4830  DATA  DELTA/0.044/ 

4840  T=0 1 

48S0  NA*Q2 

4860  MO* 03 

4870C 

4880C  CALCULATE  EG.NV.AND  EA  DIVIDED  BY  KT 
4890C 

4900  EG*< < (-3.80977E-134T  +9.95402E-10 )*T 

4910  — 8.701100E-7)*T  4. 00 003237 4l> A T  *1.155556 

4920  KT=K*T 

4930  E AP=< 0.043 8-3. 037E-8*NA»* 0.3 333 )/KT 

4940  NV*4«  82907E  1  54MD4*  1  •  54  T*4  1.5 

4950  A 00*0 .0 

496  OC 

4970C  ITERATION  TO  FIND  A  VALUE  FDR  EF  SO  THAT  P=NA- 
4980C  IT  FINOS  EF  FOR  IP-NA-)  <  0.000  1 
499 OC 
5000 
5010 
5020 
5030 
5040 
5050 
50  55 
5060 
5070 
5080 
509  0 
5100412 
5110 
5120413 


EF=0.43 

DO  12  LI  *1 .1000 

ETA=<EF-EG)/KT 

P=NV*EXP<ETA> 

IFIETA  .GE.  1. 1 P=NV4. 752254 < CETA**2 *1  .7  >**.75 > 
NAl*NA/( 1.*<4.*2.*EXP(-DELTA/KT> ) )*EXP(ETA 
fc+EAPI )+A00 
TEST=ABS(P-NAI  )  . 

DIF=P/2. 

IF<  TEST  .LE.  OIF)  GO  TO  14 
EF=EF* • 001 
CONTINUE 
WRITEC6.  13) 

F0RMATC//5X. **E  FELL  THROUGH  LOOP  0*/) 


-  i  t» 
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5130*14  EF=EF-. 001 ♦  . 00005 

5140  00  15  LI=1 .400 

5150  ETA= CEF-EG)/*T 

5160  P=NV*EXPC£TA) 

5170  IFCETA  .G£.  1.  1  P=N  V*  .  752254  I  C  E  TA*  *2  4-1  •  7  >4  *  .75  ) 

51  80  NAI=NA/  (  1.4 1  4.*2.*EXP(-DELTA/<T  )  )  )*EXPCETA* 

5185  Cr API ) ♦ADD 

5190  TEST=A8SCP-NA1 1 

5200  0  IF  — P/l 0. 

5210  I F( TE ST  .LE •  DIF)  GO  TO  30 

5220  £F=EF«-.0001 

523  0*15  CONTINUE 

5240  HP  I TE (6. 20 ) 

5250*20  FORMAT  C//5X»*HE  FELL  THROUGH  LOOP  1*/) 

5260*30  EF=EF— • 000! +.00001 

5270  00  40  L 1= 1. 400 

5280  ET  A= CEF-EG ) /KT 

5290  »=NV*EXP(ETA> 

5300  IFCETA  .GE*  l. 1P=NV*.75225*C C ETA**2* 1. 7)**. 75) 

5310  NAI=NA/C1.M4.*2.*EXP(-DELTA/KT  )  )  >*EXPCETA  + 

5315  6EAP))*AD0 

5320  TEST=ABS CP-NAII 

5330  DIF*P/100.0 

5340  I FC  TEST  .LE.  DIF1G0  TO  60 

5350  EF=EF*0. 00001 

5360*40  CONTINUE 

5370  HRITEC6.50) 

5380*50  FORMAT! //5X  »• HE  FELL  THROUGH  LOOP  2*/) 

5390*60  EF»EF— . 00001*. 000001 

5400  OO  70  Ll =1.500 

5410  5 TA=(EF-EG> /KT 

5420  P=NV*EXPCETAI 

5430  IFCETA  .G£.  1 . >P=NV* .75225* ( C ET A**2 M . 7 )**. 75) 

5440  NAI*NA/C  l.*C  4. *2. *EXP( -DELTA/XT)  >  14EXPIETA+ 

5445  tE  AP ) I +ADO 

5450  TEST=A8S(P-NAI ) 

5460  D IF=P/  1000. 

5470  IFCTEST  .LE.  DIF  ) GO  TO  90 

5480  EF=EF*0. 000001 

5490*70  CONTINUE 

5500  HRITE  C6.80) 

5510*80  FORMAT C //5X  * • HE  FELL  THROUGH  LOOP  3»/> 

5520*90  ET  A= (EF-EG )/XT 

553  OC 
5540 


PP=NAI ♦NAI *C 1.-NAI/NA) 
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55S0 

NN=ABS (NA-NAIfAOO) 

5560 

P I 0N= ( NAI —ADD! /NA 

5570 

Q4=NAI 

5580 

Q5=PP 

5590 

Q  6=NN 

5600 

Q7=P  ION 

5610 

Q8=ETA 

562  0C 

5630 

RETURN 

564  0 

END 

I 


_\r 
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S6S0C 
5710 
5660C 
5670C 
5680C 
S690C 
S69SC 
570  0C 
5720 
S730C 
574  0 
5750 
5760 
5770 
5780 
5790 
5800 
5810 
582  0 
5830 
584  0 
S8S0 
S860C 
5870 
5880 
5884 
5890 
5900 
5910 
5920 
5930C 
594  0 
5945 
5950 
596  0 
5970 
5980 
5990 
5995 
6000C 
6010 
6020 
6030 
6040 
6050 


SUBROUTINE  1NTBN0  (  Y  1 .  V2.  C  1  I.  C  22.C  12.C2I  .Cl  1 2) 

S  UBROUT I NE  TD  CALCULATE  INTERBAND  RELAXATION  TIME 
PARAMETERS  FOR  ACOUSTICAL  PHONON  SCATTERING* 

BOTH  LONGITUDINAL  AND  TRANSVERSE  VIBRATIONS 
ARE  CONSIDERED 

REAL  L21 1*L111 *L222 *L1 22  *L212 

S  TA  =  0* 2 1 25 
E  T  AM=-ET  A 
ETAS*ETA**2 
S  TAM  S=E  TAM*  *  2 
G 1=  SORT ( Y2/Y 1 ) 

G2=G1 4*2 • 

G3*G1**3. 

G4=  G1 4*4 • 

A  — ALOG  ((l.+Gl)/(l. — Gl ) ) 

9P*(  l.*G21 
BPS* { 1 .  +  G2 ) **2 • 

0MS»C l.-G2»**2. 

Pl={3./(8»*  G2 ) J  * ( 1 • +G2  —  f  BMS  *A)/( 2**G1 ) I 
s 2~ ( 3.  /( 6**G4l }*( ( 3.*G4-2.*G2+3.)/3. 
t- CBP4DMS *A)/(2.*G1 ) ) 

P3=(3./<4.*G2>  )*(8P*  ( 15.*G4-22.*G2+15.  )/(  48.4G4  1- 
1( 1  BPS*( 5.*G4- 14.*G2*5.  I /( 16.*G4> )*A/(2.*G1 ) ) » 
P4=-BMS* (1. -BP* A/(2.*G1 ) »/(4.*G2  ) 

P5=-BMS*C3.*BP-<3.+2  .*G2*3.*G4  )*A/(2.*Gl )  )/(4.*G4) 

L21 1=1 • -ETAS— 2  **ET  A*  C 1 .-ETA14P1 
&>ETA*G2* C 2.*( 1 .*ETA)*P1-3.*ETA*P2| 

L  1 1 1* ( 1 «  *ET A  1**2  *3.*ETAS/4. 

L222=(1.+FTAM)**2. *3.*ETAMS/4. 

L 122*1 .-ETA MS- 2* *E TAM* C  t.-ETAM) API* 

1 5TAM4G2*  C  2.*(1 . *ET AM ) *P 1 -3 . *ET AM *P2 ) 

L212=— ETA* ( 1 .—ETA) *P2*  G2*ET  A*  C (1  .*ETA)*P2 
t-3.*ETA*P3) 

Til  1=0.75 
T  222*0 • 75 
T 122*9. *  <  l  • — P4 1  /4 . 

T21I*T122 
T212*— 9.4P5/ 16. 
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C11*L222*G3*L122«-2.09*ET*  S*<  T2  22*G3*T  122) 
C22=U2l 1 ♦G3*L1 It *2 .09 *FT AS* ( T 2 1 l 
6*G3*TU  i)*<  Y1  **  1  •  5/Y2**  US) 

C12s|_2l2+2  • 09*ETAS*T  2 1  2  *(  Y  1**1.  5/V2**  I.  5) 
C  21  =1-2 1 2 *2. 09* ETAS* T2 1  2 

P7=3.*Gl*P5/< l.-G2)**2. 

Cl 12=P7 

RETURN 

ENO 


6060C 

60  70 
60  80 
6085 
6090 

61  00 
61  IOC 
6120 
6130 
6140C 
6150 
6160 
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